[gdal-dev] Re: C#: Retrieve boundary of coordinate system
Tomas R
monshi at home.se
Thu Jul 8 02:53:52 EDT 2010
You mean
area_of_use_code="2847"
which gives
<Area area_code="2847" area_name="Sweden - 2.5 gon W"
area_of_use="Sweden. For medium and small scale applications:- all
country. For large scale applications:- communes between approximately
14 deg 40 min and 16 deg 55 min East. See information source for map."
information_source="Lantmateriverket (National Land Survey of Sweden).
http://www.lm.se/geodesi/refsys/rt/rt_projections.html"
data_source="EPSG" revision_date="2003-06-27" deprecated="0"/>
Or is it from this line
<3021> +proj=tmerc +lat_0=0 +lon_0=15.80827777777778 +k=1.000000
+x_0=1500000 +y_0=0 +ellps=bessel +units=m +no_defs <>
or from the WKT
PROJCS["RT90 2.5 gon
V",GEOGCS["RT90",DATUM["Rikets_koordinatsystem_1990",SPHEROID["Bessel
1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],TOWGS84[414.1055246174,41.3265500042,603.0582474221,-0.8551163377,2.1413174055,-7.0227298286,0],AUTHORITY["EPSG","6124"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4124"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",15.80827777777778],PARAMETER["scale_factor",1],PARAMETER["false_easting",1500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","3021"]]
or some other file/source?
All data regarding the Swedish RT90 2.5 V where only the WKT is
complete with toWGS84 parameters.
A not so easy task, and not reliable, to retrieve the info from there I
would guess.
My goal is not to be 100% accurate at all points in the world, my goal
is to be consistent in the answers to the calls I get. That is the
functions PixelToGPS and GPSToPixel should return "compatible" results
at all parts in the world so that the main application won't have any
risk of getting confused and run into infinite loops or alike.
not to hard to step from (in WGS84) 180 to -180 and +85 to -85 and check
for every fifth degree, or so, if the transformation is valid. When
transformation not valid the plan is to switch to WGS84 which seems to
work even though results can be interesting in the boundary between the
coordinate systems.
But still, if the info already exists somewhere why not use it but I do
not know if that info exist in the data I have available.
Thanks anyway
Tomas
Francis Markham skrev 2010-07-08 08:14:
> I think the EPSG dictionary defines areas of validity for the defined
> coordinate systems. That could be a starting point, although there
> may be a better solution
>
> -Francis
>
> On 8 July 2010 16:01, Tomas R<monshi at home.se> wrote:
>> No one?
>>
>> Ok, I guess then that there are no common function/method of doing this. My
>> approach will the be trial and error. From WGS 84 transform to the local
>> coordinate system and back again to WGS84 and check if the result seems
>> valid.
>>
>> Should work or any expert here have another solution?
>>
>> Yours
>> Tomas
>>
>> Tomas R skrev 2010-07-05 09:05:
>>> Attacking GDAL via the C# interface and I have a question
>>>
>>> Have set up two (Osr) SpatialReferences and transformations to and form
>>> both. One system is WGS84 and the other is unknown.
>>>
>>> To get a consistent behaviour I need to retrieve the boundary, in WGS84
>>> coordinates, of the unknown system. If known I can fallback on WGS84 when
>>> calls outside the scope of the coordinates system are made.
>>> Should be easy to read of some attribute of the SpatialReference or
>>> alike. But which/how?
>>>
>>> Any suggestions?
>>>
>>> Or how should I solve it - in general words:
>>> My application maps coordinates from and to WGS84 and serves them to
>>> another via calls like PixelToGPS and GPSToPixel. The originating coordinate
>>> system is always WGS84. Need to handle calls made outside the scope of the
>>> coordinate system and my approach, to get a consisting behaviour of the two
>>> functions, is to revert to WGS84 which is defined at all possible points.
>>>
>>>
>>> Yours
>>> Tomas
>> _______________________________________________
>> gdal-dev mailing list
>> gdal-dev at lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>>
More information about the gdal-dev
mailing list