[Proj] Problem overlaying georeferenced images in Google, Maps(projection problem!)
Jorge
jorge.arevalo at gmail.com
Fri May 30 04:12:20 PDT 2008
Now, the results in Google Earth:
The image before reproyecting (UTM zone 30 datum ED50)
http://www.nebulared.com/tmp_geo/Balearic_Island_UTM_ED50.png
The image reproyected
http://www.nebulared.com/tmp_geo/Balearic_Island_Reprojected.png
What's the problem?:
- The input image's projection really isn't UTM ED50 (Mapinfo lies??), or is bad georeferenced (low probability...)
- The output image's projection really isn't GMaps projection (assuming GMaps projection = GEarth projection)
- GMaps projection not equal GEarth projection (I think so...)
- Any other?
The line executed:
gdalwarp
-s_srs "+init=epsg:23030 +towgs84=-84,-107,-120,0,0,0,0"
-t_srs "+init=epsg:3785 +towgs84=0,0,0,0,0,0,0"
-tr 250 250
-co INTERLEAVE=PIXEL
Balearic_Island_UTM_ED50.png
Balearic_Island_Reprojected.png
The epsg codes:
Input code (23030):
# ED50 / UTM zone 30N
<23030> +proj=utm +zone=30 +ellps=intl +units=m +no_defs <>
Output code (3785):
# Spherical/Web Mercator
<3785> +proj=merc +lat_ts=0.0 +lon_0=0.0 +k=1.0 +x_0=0.0 +y_0=0.0 +a=6378137.0 +b=6378137.0 +units=m +nadgrids=@null +no_defs <>
Sorry for the insistence, and thanks a lot
Jorge
Jorge wrote:
> The EPSG code 3785 is a custom one, but the same of 900913 (I didn't know it thanks), except for the "+wktext" term. What does it mean?
>
> OK, I'm lost with this issue. Could my calcs be wrong? I have, for example, this image:
> http://www.nebulared.com/tmp_geo/UTM_ED50_Balearic_Island.png
>
> And is georeferenced by this world file:
> http://www.nebulared.com/tmp_geo/UTM_ED50_Balearic_Island.wld
>
> The image is obtained for a MIF/MID couple of files with UTM projection, zone 30, ED50 datum (Mapinfo said... I don't know what
> are the dX, dY, dZ parameters, buth they must be -84, -107, -120). Really, the raster file obtained from MIF/MID files isn't a
> PNG file, but I did the transform to "see" what I was doing.
>
> OK, I need to reproject the image, and the output image must be used to overlay it in a GMaps' frame. This is my last attempt:
>
> gdalwarp
> -s_srs "+init=epsg:23030 +towgs84=-84,-107,-120,0,0,0,0"
> -t_srs "+init=epsg:3785 +towgs84=0,0,0,0,0,0,0"
> -tr 250 250
> -co INTERLEAVE=PIXEL
> UTM_ED50_Balearic_Island.png
> GMaps_Balearic_Island.tif
>
> Then, I transform the TIFF file in a PNG (transformation OK) and I have this image:
> http://www.nebulared.com/tmp_geo/GMaps_Balearic_Island.png
>
> But the image doesn't fit exactly in the map... Is there any error in this reprojection? If the projection is OK, I need to review
> the rest of the code...
More information about the Proj
mailing list