[gdal-dev] write error
Even Rouault
even.rouault at mines-paris.org
Sun Apr 11 05:19:38 EDT 2010
The RasterIO parameters are not correct.
Try instead :
poSrcBand->RasterIO(GF_Read,0,i,nX,1,pafScanline,nX,1,GDT_Float32,0,0);
poDstBand->RasterIO(GF_Write,0,i,nX,1,pafScanline,nX,1,GDT_Float32,0,0);
The 'for(int i=1;i<nY;i++)' loop is not correct. The starting value of i
should be 0.
Ther's also a leak of the *pafScanline buffer.
But if I were you, I'd use the following function to do the copy of the
imagery :
CPLErr CPL_DLL CPL_STDCALL GDALDatasetCopyWholeRaster(
GDALDatasetH hSrcDS, GDALDatasetH hDstDS, char **papszOptions,
GDALProgressFunc pfnProgress, void *pProgressData );
Le Sunday 11 April 2010 08:32:51 weixj2003ld, vous avez écrit :
> I use the following code to read data from a geotiff file and write it into
> a new one.the bands number is 6; ...
> GDALDataset *poSrcDS,*poDstDS;//poSrcDs:the source file,poDstDS:the goal
> source GDALRasterBand *poSrcBand,*poDstBand;//poSrcBand:the source
> band,poDstBand:the goal band GDALDriver *poDstDriver;
> double *m_AdGeoTransform=new double[6];
>
> ...
> poDstDriver=GetGDALDriverManager()->GetDriverByName(poSrcDS->GetDriver()->G
>etDescrition());
> poDstDS=poDstDriver->create(DstFilename,nX,nY,6,GDT_Float32,NULL);
> poSrcDS->GetGeoTransform(m_AdGeoTransform)
> poDstDS->SetGeoTransform(m_AdGeoTransform);
> poDstDS->SetProjection(poSrcDS->GetProjectionRef());
>
> for(int iBand=0;iBand<6;iBand++)
> {
> poSrcBand=poSrcDS->GetRaserBand(iBand+1);
> poDstBand=poDstDS->GetRaserBand(iBand+1);
> float *pafScanline;
> pafScanline=(float *)CPLMalloc(sizeof(float)*nX;
> for(int i=1;i<nY;i++)
> {
>
> poSrcBand->RasterIO(GF_Read,0,0,nX,i,pafScanline,nX,1,GDT_Float32,0,0);
> poDstBand->RasterIO(GF_Write,0,0,nX,i,pafScanline,nX,1,GDT_Float32,0,0);//I
>f I comment this code,the program will be ok. }
> }
> GDALClose((GDALDatasetH)poDstDS);
>
> When I run my program,CPU is very alway busy,and the program is
> blocked,why?
More information about the gdal-dev
mailing list