Azzera filtri
Azzera filtri

Validation and Correction of GPS Altitude with data out of Digital Elevation model

5 visualizzazioni (ultimi 30 giorni)
Hello.
I have a Project which also consist of validating and correcting GPS Altitude data.
I have a 2Hz GPS log and i want to verify the Altitude through Latitude and Longitude in a DEM.
The 10m resolution DEM is from a government webpage and is an ASCII-File (coordinate system ETRS89/WGS84).
(https://www.data.gv.at/katalog/dataset/land-stmk_digitalesgelndemodell10m/resource/ed0d0239-8c06-4e51-b1f1-59001e23dcc9
So the data is present in meter x- and y-coordinates not in degree lat and lon.
I successfully read the data using:
[A,R] = readgeoraster('ALS_DGM_10M_UTM33N.asc','OutputType','double');
That outputs R as:
R =
MapCellsReference with properties:
XWorldLimits: [388131.5 589511.5]
YWorldLimits: [5161133.5 5298053.5]
RasterSize: [13692 20138]
RasterInterpretation: 'cells'
ColumnsStartFrom: 'north'
RowsStartFrom: 'west'
CellExtentInWorldX: 10
CellExtentInWorldY: 10
RasterExtentInWorldX: 201380
RasterExtentInWorldY: 136920
XIntrinsicLimits: [0.5 20138.5]
YIntrinsicLimits: [0.5 13692.5]
TransformationType: 'rectilinear'
CoordinateSystemType: 'planar'
I found the program "deg2utm" which converts my polar coordinates in to cartesian coordinates. (which it does successfully)
e.g.
[x_GPS,y_GPS,utmzone] = deg2utm(lat_GPS,lon_GPS)
However neither geointerp nor mapprofile works with my "R" which is "map.rasterref.MapCellsReference"
Error using geointerp
Expected R to be one of these types:
map.rasterref.GeographicCellsReference, map.rasterref.GeographicPostingsReference
Instead its type was map.rasterref.MapCellsReference.
Error in geointerp (line 22)
validateattributes(R, ...
Error using internal.map.convertToGeoRasterRef (line 39)
Function MAPPROFILE expected input argument 2, 'R', to be a 3-element referencing vector, a 3-by-2 referencing matrix, or a (scalar) geographic raster
reference object.
Error in mapprofile (line 123)
R = internal.map.convertToGeoRasterRef( ...
Maybe somebody can help me resolve my issues
THANKS

Risposte (0)

Prodotti


Release

R2020a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by