[ gpsdrive ] Map projection

Ryan J. Earley ryan.j.earley at gmail.com
Sat Aug 4 08:02:17 AKDT 2007


On 8/2/07, Petra Galuscakova <galuscakova at gmail.com> wrote:
> Hello,
> I would like to ask you, in which projection  the maps that are
> displayed in Gpsdrive are. I`ve found, that it`s Miller projection.
> Could you confirm that?
>
> I would like to use gdalwarp software to translate my map into the
> projection which Gpsdrive maps uses, but I`m not sure how to do it. My
> map is in S-JTSK projection (ESRI 102067).
>  I have tried to translate it into Miller projection (EPSG 53002).
> To do that I have used next:  gdalwarp -s_srs "+init=esri:102067"
> -t_srs "+init=epsg:53002" map_file0025g.gif map_file0025t.gif
> I`ve compared this map with expedia map from the same area and the
> same scale (downloaded by Gpsdrive) and I found some differences.
> Am I doing something wrong or are there possible some differences
> between maps?  Am I using appropriate epsg number?

As long as gdal understands your projection, you can use gdalwarp to
go from your projection to WGS84/LAT-LON (EPSG:4326).  Then you should
be able to treat the resulting maps as "top_" files in gpsdrive.

Example: i recently had a dataset in NJ State Plane i needed in
gpsdrive.  I used the following with  success:
gdalwarp -s_srs "+init=nad83:2900 +datum=NAD83 +units=us-ft \
+proj=tmerc +lon_0=-74d30 +lat_0=38d50 +k=.9999 +x_0=150000 +y_0=0 \
+no_defs +ellps=GRS80 +towgs84=0,0,0" -t_srs "+init=epsg:4326 \
+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0" \
INFILE.geotif     OUTFILE.latlon.gtiff

if you replace the -s_srs parameters (inside the quotes) with your
map's parameters, this conversion should "just work" and will give you
a gpsdrive "top" file.

gpsdrive has 2 (so far?) projections -- "map" and "top".  here is one
explanation:

<Cut-and-paste from gdal_slice.sh
http://hamish.bowman.googlepages.com/gdal_slice.sh >
# Important! The maps must be named "map_*" for UTM-like projections
# (lat:lon = 1:cos(lat)) and "top_*" for lat/lon Plate carree projection
# (lat:lon = 1:1). The prefix is given so that gpsdrive knows how to
# scale the maps correctly. Alternatively the maps can be stored without
# prefix in subdirectories of $HOME/.gpsdrive/ which end in "_map" or
# "_top".
</Cut-and-paste>

If you original map is not 1280x1024px, you can use Hamish's
gdal_slice.sh, a geotiff-->gpsdrive tiler, to generate the gpsdrive
tiles and map_koord entries quickly.

it is available here:
http://hamish.bowman.googlepages.com/gpsdrivefiles

hope that helps,
-ryan


More information about the GPSdrive mailing list