Help with rasters

[Jeff Dege]

I've been trying to do the same, but the shp2img that came with ms4w
won't output TIFF files.

> -----Original Message-----
> From: Brent Pedersen [mailto:bpederse at nature.Berkeley.EDU] 
> Sent: Friday, October 20, 2006 12:38 PM
> To: Jeff Dege
> Cc: MAPSERVER-USERS at LISTS.UMN.EDU
> Subject: Re: [UMN_MAPSERVER-USERS] Help with rasters
> 
> hi, i just did something similar with a bunch of MrSID imagery,
> i got a bunch of help from #gdal,
> 
> i generated a 15K * 15K pixel (.tiff) overview using shp2img, 
> then (after gdaladdo) used that for large scales in the .map file
> and only resort to the MrSIDs when zoomed in pretty far.
> 
> by the way, i think it took > 1 hour to generate the 15K by 15K
> pixel overview. but once it's done...
> 
> 
> -brent
> 
> 
> 
> On Fri, 20 Oct 2006, Jeff Dege wrote:
> 
> > Yep. That did it. I know have an image of the full extent 
> of the map.
> >
> > Only problem is that it took ten minutes to generate.
> >
> > Next?  Read up on optimization and tuning of raster maps.
> >
> > Thanks for the help.
> >
> >
> >> -----Original Message-----
> >> From: Bart van den Eijnden (OSGIS) [mailto:bartvde at osgis.nl]
> >> Sent: Friday, October 20, 2006 11:49 AM
> >> To: Jeff Dege
> >> Cc: MAPSERVER-USERS at LISTS.UMN.EDU
> >> Subject: Re: [UMN_MAPSERVER-USERS] Help with rasters
> >>
> >> First of all, change:
> >>
> >> DATA "rasterindex"
> >>
> >> into:
> >>
> >> TILEINDEX "rasterindex"
> >>
> >> Best regards,
> >> Bart
> >>
> >> Jeff Dege schreef:
> >>> A customer has sent us a bunch of mapdata, to see if we 
> can make any
> >>> sense from it.  I've been messing around with mapserver 
> for a while,
> >>> now, and I've figured out how to work some of its
> >> capabilities, but this
> >>> is my first attempt to handle rasters.
> >>>
> >>> The data set has four sets of shape files, three of which are
> >>> point/line/polygon data and the fourth which seems to be a
> >> tileindex for
> >>> the thousands of jpg and jgw files.
> >>>
> >>> I'm starting simple - writing a map file and trying to get
> >> a map image
> >>> using shp2img.  I have no problem doing this with the three
> >>> point/line/polygon layers.  But I've not been able to do
> >> this with the
> >>> raster data.
> >>>
> >>> I can run gdalinfo on the jpg files:
> >>>
> >>>   Driver: JPEG/JPEG JFIF
> >>>   Size is 1610, 1609
> >>>   Coordinate System is `'
> >>>   Origin = (-93.112361,44.933805)
> >>>   Pixel Size = (0.00001279,-0.00000900)
> >>>   Corner Coordinates:
> >>>   Upper Left  ( -93.1123608,  44.9338046)
> >>>   Lower Left  ( -93.1123608,  44.9193245)
> >>>   Upper Right ( -93.0917678,  44.9338046)
> >>>   Lower Right ( -93.0917678,  44.9193245)
> >>>   Center      ( -93.1020643,  44.9265646)
> >>>   Band 1 Block=1610x1 Type=Byte, ColorInterp=Red
> >>>   Band 2 Block=1610x1 Type=Byte, ColorInterp=Green
> >>>   Band 3 Block=1610x1 Type=Byte, ColorInterp=Blue
> >>>
> >>> I can run gdaltindex to create my own tileindex shape file:
> >>>
> >>>   gdaltindex  \\server\share\test\data\rasterindex
> >>> \\server\share\test\data\data\*.jpg
> >>>
> >>> And I can run ogrinfo on the generated shapefile:
> >>>
> >>>   INFO: Open of `\\server\share\app\data'
> >>>   using driver `ESRI Shapefile' successful.
> >>>
> >>>   Layer name: rasterindex
> >>>   Geometry: Polygon
> >>>   Feature Count: 683
> >>>   Extent: (-93.338884, 44.455961) - (-92.721094, 44.933805)
> >>>   Layer SRS WKT:
> >>>   (unknown)
> >>>   location: String (255.0)
> >>>
> >>> But when I run shp2img, I get nada.  With DEBUG=ON:
> >>>
> >>>   Unable to open file rasterindex for layer A ... ignoring
> >> this missing
> >>> data.
> >>>   msDrawMap(): Layer 0 (A), 0.000s
> >>>   msDrawMap(): Drawing Label Cache, 0.000s
> >>>   msDrawMap() total time: 0.000s
> >>>
> >>> The rastedindex shape files (.dbf, .shp, and .shf) files
> >> are in the same
> >>> directories as the other point/line/polygon shape files
> >> that I've been
> >>> able to successfully draw layers with using this mapfile.
> >> Ogrinfo has
> >>> no trouble accessing the file.
> >>>
> >>> The underlying .jpg files are in the same directory, and
> >> gdalinfo and
> >>> gdaltindex had no trouble opening them.
> >>>
> >>> The layer in the mapfile itself is pretty simple:
> >>>
> >>>   LAYER
> >>>     NAME "A"
> >>>     TYPE RASTER
> >>>     STATUS ON
> >>>     DATA "rasterindex"
> >>>   END
> >>>
> >>> Any hints as to what I'm doing wrong?
> >>>
> >>>
> >>>
> >>
> >>
> >> --
> >> Bart van den Eijnden
> >> OSGIS, Open Source GIS
> >> bartvde at osgis.nl
> >> http://www.osgis.nl
> >>
> >>
> >
>