[mapserver-users] Projection issue with vrt file

Fri Apr 20 21:39:17 EDT 2012

Hi All,

I'm trying to add a VRT raster file (BSB chart) to my mapfile, but it doesn't work totally.
In fact, all of my tiles are generated, but they have no data...  No error is logged...

It sounds like the projection or the transformation matrix is not well supported, or not well configured.

If I put in comments the GeoTransform element, some blank tiles are shown where the lower left corner is located at the origin (0.000, 0.000).

Moreover, I successfully opened that same VRT file with QGIS, and the projection is ok...

A copy of my VRT file is shown just after my message...

Here is a snippet of my mapfile:

  LAYER
    NAME "131701"
	GROUP "BSB"
	#EXTENT -70.980392 46.777300 -70.191770 47.279794
    	DATA "BSB/131701.vrt"
	TYPE RASTER
    	STATUS OFF
	CLASS
      		NAME 'BSB'
    	END
    	METADATA
      		WMS_TITLE "BSB"
	  WMS_SRS "EPSG:4326 EPSG:29450 EPSG:2949 EPSG:32187 EPSG:32188 EPSG:900913"
    	END	
  END

Thanks in advance for your support,

Marc-Andre

---------------------

<VRTDataset rasterXSize="11388" rasterYSize="8091">
  <SRS>+proj=merc +lat_ts=47.166667 +datum=NAD83 +nodefs</SRS>
<GeoTransform>-5.3802490306424974e+006, 5.0001511315648770e+000,-1.8211665160056521e-004, 4.0492551038299361e+006,-2.5595698020313497e-004,-5.0010081232801378e+000</GeoTransform>
  <Metadata>
    <MDI key="BSB_KNP">SC=50000,GD=NAD83,PR=MERCATOR,PP=47.166667,PI=5.000000,SP=,SK=0.000000,TA=90.000000,UN=METRES,SD=LNT,DX=5.000000,DY=5.000000</MDI>
    <MDI key="BSB_KNQ">EC=RF,GD=NAR,VC=ZZZ,SC=ZZZ,PC=MC,P1=-70.582029,P2=47.166667,P3=0.0,P4=NOT_APPLICABLE,GC=GE,RM=NONE</MDI>
    <MDI key="CUTLINE">MULTIPOLYGON(((257.76284645559468 7745.7076975815671,2490.025011025306 7746.5301512821934,2492.2163879532559 7181.6425479723066,3276.1412438760321 7210.5951029770404,4429.0270760545409 6817.9091369725729,4890.9737230091405 6736.5650032022368,11250.811878184026 354.19735628215676,11251.785400765162 244.45756921568665,3537.7418346315158 244.33482436588929,3523.5069250460119 2029.8367618071302,1881.7430005409099 2028.8818991008357,1868.6640235103462 2447.8575431625686,258.92730402280057 2434.9666589315798)))</MDI>
    <MDI key="DESCRIPTION">SAULT-AU-COCHON À/TO QUÉBEC</MDI>
  </Metadata>
  <GCPList Projection="+proj=merc +lat_ts=47.166667 +datum=NAD83 +nodefs">
    <GCP Id="1" Pixel="259.0000" Line="7746.0000" X="-5.378954975641E+006" Y="4.010515470308E+006" />
    <GCP Id="2" Pixel="4892.0000" Line="7746.0000" X="-5.355789271702E+006" Y="4.010515470308E+006" />
    <GCP Id="3" Pixel="4892.0000" Line="6736.0000" X="-5.355789271702E+006" Y="4.015564569030E+006" />
    <GCP Id="4" Pixel="11252.0000" Line="354.0000" X="-5.323986290417E+006" Y="4.047476941272E+006" />
    <GCP Id="5" Pixel="11252.0000" Line="244.0000" X="-5.323986290417E+006" Y="4.048033445911E+006" />
    <GCP Id="6" Pixel="259.0000" Line="244.0000" X="-5.378954975641E+006" Y="4.048033445911E+006" />
    <GCP Id="7" Pixel="9527.0000" Line="986.0000" X="-5.332614722193E+006" Y="4.044323691092E+006" />
    <GCP Id="8" Pixel="7000.0000" Line="986.0000" X="-5.345251249970E+006" Y="4.044323691092E+006" />
    <GCP Id="9" Pixel="4472.0000" Line="986.0000" X="-5.357887777748E+006" Y="4.044323691092E+006" />
    <GCP Id="10" Pixel="1945.0000" Line="986.0000" X="-5.370524305526E+006" Y="4.044323691092E+006" />
    <GCP Id="11" Pixel="8263.0000" Line="2837.0000" X="-5.338932986081E+006" Y="4.035066484912E+006" />
    <GCP Id="12" Pixel="5736.0000" Line="2837.0000" X="-5.351569513859E+006" Y="4.035066484912E+006" />
    <GCP Id="13" Pixel="3208.0000" Line="2837.0000" X="-5.364206041637E+006" Y="4.035066484912E+006" />
    <GCP Id="14" Pixel="681.0000" Line="2837.0000" X="-5.376842569415E+006" Y="4.035066484912E+006" />
    <GCP Id="15" Pixel="4472.0000" Line="4685.0000" X="-5.357887777748E+006" Y="4.025823838111E+006" />
    <GCP Id="16" Pixel="1945.0000" Line="4685.0000" X="-5.370524305526E+006" Y="4.025823838111E+006" />
    <GCP Id="17" Pixel="3209.0000" Line="6531.0000" X="-5.364206041637E+006" Y="4.016595685774E+006" />
    <GCP Id="18" Pixel="682.0000" Line="6531.0000" X="-5.376842569415E+006" Y="4.016595685774E+006" />
  </GCPList>
  <VRTRasterBand dataType="Byte" band="1">
    <Metadata />
    <ColorInterp>Palette</ColorInterp>
    <ColorTable>
      <Entry c1="0" c2="0" c3="0" c4="255" />
      <Entry c1="255" c2="255" c3="255" c4="255" />
      <Entry c1="255" c2="214" c3="140" c4="255" />
      <Entry c1="107" c2="165" c3="222" c4="255" />
      <Entry c1="182" c2="213" c3="240" c4="255" />
      <Entry c1="100" c2="160" c3="140" c4="255" />
      <Entry c1="0" c2="86" c3="158" c4="255" />
      <Entry c1="75" c2="75" c3="75" c4="255" />
      <Entry c1="123" c2="8" c3="90" c4="255" />
      <Entry c1="100" c2="100" c3="100" c4="255" />
      <Entry c1="243" c2="137" c3="71" c4="255" />
      <Entry c1="44" c2="99" c3="207" c4="255" />
      <Entry c1="245" c2="255" c3="80" c4="255" />
      <Entry c1="255" c2="0" c3="0" c4="255" />
    </ColorTable>
    <SimpleSource>
      <SourceFilename relativeToVRT="0">E:\ms4w\apps\marinfo\data\BSB\131701.kap</SourceFilename>
      <SourceBand>1</SourceBand>
      <SourceProperties RasterXSize="11388" RasterYSize="8091" DataType="Byte" BlockXSize="11388" BlockYSize="1" />
      <SrcRect xOff="0" yOff="0" xSize="11388" ySize="8091" />
      <DstRect xOff="0" yOff="0" xSize="11388" ySize="8091" />
    </SimpleSource>
  </VRTRasterBand>
</VRTDataset>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/mapserver-users/attachments/20120420/663138b5/attachment.html