Hej,
Jag har ett litet problem som jag hoppas att någon kunnig själ kan hjälpa mig med!
Jag håller på och skapa en webb-applikation som ska visualisera gps-tracks på rasterkartor, mha mapserver.
Jag har inga problem när det gäller att få RIK-kartor (Gröna kartan tex) att fungera, men värre är det med BSB-filer (Sjökort - Ostkusten).
Jag får inte riktigt till det med att definiera rätt projektioner och datum.
Om man ser till papperssjökort så är dom ritade med Gauss (Transversal Mercator) projektion och WGS84 datum. Det är ju iaf vad som står på själva korten.
Och om jag har förstått det rätt så är det projektionen RT90-2,5gonV som avses här.
Om jag däremot kontrollerar den information som finns i .KAP filerna, så ska dom vara definierade i följande proj/datum:
Driver: BSB/Maptech BSB Nautical Charts
Size is 13347, 9213
Coordinate System is `'
GCP Projection = GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG",7030]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG",6326]],PRIMEM["Greenwich",0,AUTHORITY["EPSG",8901]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG",9108]],AXIS["Lat",NORTH],AXIS["Long",EAST],AUTHORITY["EPSG",4326]]
Med andra ord longlat WGS84.
Och det verkar stämma bra, om man ritar ut dom utan att projicera om, så får man en bild som är lite för bred, men i övrigt korrekt. Vilket känns logiskt när man ritar ut lika många grader på båda axlarna.
Om jag istället försöker att projicera om det hela, och rita ut i RT90, så kan jag inte få det att stämma med medel-meridianen. Det resulterande kortet lutar alltid lite åt något håll, även om jag sätter medelmeridianen till kortets egna medelmeridian.
Så långt mitt problem.
Jag antar att det finns en uppsättning projektionsparametrar & datum som stämmer, eller så har jag helt enkelt gjort någon annan tanke-vurpa.
Det konstiga är att allt fungerar så bra med RIK-filerna bara.
Uppskattar alla svar!!
Med vänlig hälsning
Marcus Kempe
-------------------
Mapserver använder PROJ.4 för transformeringar och projektioner. Här är dom jag tycker borde fungera, men som skapar ett resultat som lutar ett antal grader.
output projektion:
PROJECTION
"init=epsg:3021" # RT90 2.5 gon V
# <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 <>
END
och input:
PROJECTION
"init=epsg:4326" # WGS84
END
Jag har ett litet problem som jag hoppas att någon kunnig själ kan hjälpa mig med!
Jag håller på och skapa en webb-applikation som ska visualisera gps-tracks på rasterkartor, mha mapserver.
Jag har inga problem när det gäller att få RIK-kartor (Gröna kartan tex) att fungera, men värre är det med BSB-filer (Sjökort - Ostkusten).
Jag får inte riktigt till det med att definiera rätt projektioner och datum.
Om man ser till papperssjökort så är dom ritade med Gauss (Transversal Mercator) projektion och WGS84 datum. Det är ju iaf vad som står på själva korten.
Och om jag har förstått det rätt så är det projektionen RT90-2,5gonV som avses här.
Om jag däremot kontrollerar den information som finns i .KAP filerna, så ska dom vara definierade i följande proj/datum:
Driver: BSB/Maptech BSB Nautical Charts
Size is 13347, 9213
Coordinate System is `'
GCP Projection = GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG",7030]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG",6326]],PRIMEM["Greenwich",0,AUTHORITY["EPSG",8901]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG",9108]],AXIS["Lat",NORTH],AXIS["Long",EAST],AUTHORITY["EPSG",4326]]
Med andra ord longlat WGS84.
Och det verkar stämma bra, om man ritar ut dom utan att projicera om, så får man en bild som är lite för bred, men i övrigt korrekt. Vilket känns logiskt när man ritar ut lika många grader på båda axlarna.
Om jag istället försöker att projicera om det hela, och rita ut i RT90, så kan jag inte få det att stämma med medel-meridianen. Det resulterande kortet lutar alltid lite åt något håll, även om jag sätter medelmeridianen till kortets egna medelmeridian.
Så långt mitt problem.
Jag antar att det finns en uppsättning projektionsparametrar & datum som stämmer, eller så har jag helt enkelt gjort någon annan tanke-vurpa.
Det konstiga är att allt fungerar så bra med RIK-filerna bara.
Uppskattar alla svar!!
Med vänlig hälsning
Marcus Kempe
-------------------
Mapserver använder PROJ.4 för transformeringar och projektioner. Här är dom jag tycker borde fungera, men som skapar ett resultat som lutar ett antal grader.
output projektion:
PROJECTION
"init=epsg:3021" # RT90 2.5 gon V
# <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 <>
END
och input:
PROJECTION
"init=epsg:4326" # WGS84
END