Validation and Correction of GPS Altitude with data out of Digital Elevation model
2 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
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
0 comentarios
Respuestas (0)
Ver también
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!