[Qgis-user] projection change

Micha Silver micha at arava.co.il
Thu Dec 17 04:20:33 PST 2009


On 17/12/2009 13:27, ahmettemiz88 at gmail.com wrote:
> thank you again
>
>
> Original map is in *adf format it's metadata like this:
>
> name: Lambert Conformal Conic
> proj: lcc
> datum: wgs84
> ellps: wgs84
> lat_1: 40.66666666666666
> lat_2: 43.33333333333334
> lat_0: 0
> lon_0: 34
> x_0: 1000000
> y_0: 0
> no_defs: defined
OK, do you have other layers in that CRS? Do they align with the new 
shapefile??
>
> I save it as shp fie in QGIS.
And what CRS did you specify when you saved as shapefile (I can't find 
any CRS in QGIS like the one above...). In Settings->Options->CRS 
there's a setting how QGIS treats a layer when it cannot recognize the 
CRS. If you have it set to some global value, then QGIS will assume the 
new shapefile should be in that CRS, and set it as such. Have a look in 
the shapefile's metadata (property page). What CRS does it say there?

Also, I think that ogr2ogr will read the adf file directly (is it an 
Arc/Info Binary?). Then you won't need the interim step of converting to 
a shapefile.

-- 
Micha
>
> then I am trying to convert it lat/long with ogr2ogr as:
>
> ogr2ogr -t_srs EPSG:4326 -s_srs '+proj=lcc +lat_1=40.66666666 
> +lat_2=43.333333333 +lat_0=0 lon_0=34.0 +x_0=1000000 +y_0=0 
> +datum=wgs84 +spheroid=wgs84 +units=m +no_defs' 
> geohidro_cizgi_adf2.shp hidro_cizgi_adf.shp
>
> 17 Ara 2009 12:50 tarihinde, Micha Silver <micha at arava.co.il> şunu yazdı:
> > On 17/12/2009 10:49, ahmet temiz wrote:
> >
> >
> > thank you
> >
> >
> >
> > I tried as you have suggested. But there seems to be something wrong.
> >
> >
> >
> > resulting map's extent looks weird.
> >
> >
> >
> >
> > I agree with you on that...
> >
> > Are you sure about the original shapefile projection? Does it have 
> an attached *.prj file? If yes, then you should be able to run ogr2ogr 
> *without* specifying the -s_srs.
> >
> >
> >
> > I can't find in the EPSG database any projection with your 
> parameters. Is it some special CRS?
> >
> >
> >
> >
> > $ ogr2ogr -t_srs EPSG:4326 -s_srs  '+proj=lcc +lat_1=40.66666666
> >
> > +lat_2=43.333333333 +lat_0=41 lon_0=34.0 +x_0=1000000 +y_0=0
> >
> > +datum=wgs84 spheroid=wgs84 +units=m' geohidro_cizgi_adf2.shp
> >
> > hidro_cizgi_adf.shp
> >
> > orkun at orkun-desktop:~$
> >
> > orkun at orkun-desktop:~$
> >
> > orkun at orkun-desktop:~$ ogrinfo geohidro_cizgi_adf2.shp -al -summary
> >
> > INFO: Open of `geohidro_cizgi_adf2.shp'
> >
> >       using driver `ESRI Shapefile' successful.
> >
> >
> >
> > Layer name: geohidro_cizgi_adf2
> >
> > Geometry: Line String
> >
> > Feature Count: 4929
> >
> > Extent: (-28.672066, 75.100678) - (37.866883, 81.508631)  *****
> >
> > Layer SRS WKT:
> >
> > GEOGCS["GCS_WGS_1984",
> >
> >     DATUM["WGS_1984",
> >
> >         SPHEROID["WGS_1984",6378137,298.257223563]],
> >
> >     PRIMEM["Greenwich",0],
> >
> >     UNIT["Degree",0.017453292519943295]]
> >
> > UserId: Integer (10.0)
> >
> > FNODE_: Integer (10.0)
> >
> > TNODE_: Integer (10.0)
> >
> > LPOLY_: Integer (10.0)
> >
> > RPOLY_: Integer (10.0)
> >
> > LENGTH: Real (32.3)
> >
> > HIDRO_CIZG: Integer (10.0)
> >
> > HIDRO_CIZG: Integer (10.0)
> >
> > F_CODE: String (80.0)
> >
> > F_NAME: String (80.0)
> >
> > SYMBOL: Integer (10.0)
> >
> > DISPSCALE: Integer (10.0)
> >
> >
> >
> >
> >
> > 2009/12/16 Micha Silvermicha at arava.co.il>:
> >
> >
> >
> >
> > ahmettemiz88 at gmail.com wrote:
> >
> >
> >
> >
> >
> >
> > hello
> >
> >
> >
> > I have difficulty changing map projection from lamber conical to
> >
> > geographic ?
> >
> >
> >
> > how can I do that ?
> >
> >
> >
> > here is the original map's metada:
> >
> >
> >
> > name: Lambert Conformal Conic
> >
> > proj: lcc
> >
> > datum: wgs84
> >
> > ellps: wgs84
> >
> > lat_1: 40.66666666666666
> >
> > lat_2: 43.33333333333334
> >
> > lat_0: 0
> >
> > lon_0: 34
> >
> > x_0: 1000000
> >
> > y_0: 0
> >
> > no_defs: defined
> >
> >
> >
> > I tried with ogr2ogr :
> >
> >
> >
> > ogr2ogr -t_srs '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs' -s_srs
> >
> > '+proj=lcc +lat_1=40.66666666 +lat_2=43.333333333 +lat_0=0 lon_0=34
> >
> > +x_0=1000000 +y_0=0 +datum=wgs84 +units=meter' geohidro_cizgi_adf94.shp
> >
> > hidro_cizgi_adf.shp
> >
> >
> >
> >
> >
> >
> >
> >
> > I think you might be having trouble with the quoting. Maybe try with 
> just
> >
> > the EPSG code, i.e.:
> >
> > ogr2ogr -t_srs EPSG:4326 geohidro_cizgi_adf94.shp hidro_cizgi_adf.shp
> >
> >
> >
> > --
> >
> > Micha
> >
> >
> >
> >
> > regards
> >
> >
> >
> > Ahmet Temiz
> >
> >
> >
> > This mail was received via Mail-SeCure System.
> >
> >
> >
> >
> >
> >
> >
> > ------------------------------------------------------------------------
> >
> >
> >
> > _______________________________________________
> >
> > Qgis-user mailing list
> >
> > Qgis-user at lists.osgeo.org
> >
> > http://lists.osgeo.org/mailman/listinfo/qgis-user
> >
> >
> >
> > This mail was received via Mail-SeCure System.
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> > This mail was received via Mail-SeCure System.
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> This mail was received via Mail-SeCure System.




More information about the Qgis-user mailing list