read a Quickbird multispectral image

Question

0.00/5 (No votes)

See more:

I am trying to read a quickbird multispectral image which is having 16 bits pixel value each for red, green, blue and alpha component. I am reading its pixel values using gdal library but I am not getting the exact values of all 4 layers. I am getting very different image pixel values than the original and I dont know why it is coming so. please help me.........

Below is my code to read quickbird multispectral image....

C++

GDALDataset *hd1; 
GDALRasterBand *Rband, *Gband, *Bband, *Alphaband; 
int widthMS, 
int *Rbuffer, *Gbuffer, *Bbuffer, *Alphabuffer; 
int i,j,k; 
hd1 =(GDALDataset *)GDALOpen("G:\\images\\QB-multispectral-clipd.tif",GA_ReadOnly); 
widthMS = hd1->GetRasterXSize(); 
heightMS = hd1->GetRasterYSize(); 
Rband = hd1->GetRasterBand(1); 
Gband = hd1->GetRasterBand(2); 
Bband = hd1->GetRasterBand(3); 
Alphaband = hd1->GetRasterBand(4); 
Rbuffer = (int*)CPLMalloc(sizeof(int)*widthMS*heightMS); 
Gbuffer = (int*)CPLMalloc(sizeof(int)*widthMS*heightMS); 
Bbuffer = (int*)CPLMalloc(sizeof(int)*widthMS*heightMS); 
Alphabuffer = (int*)CPLMalloc(sizeof(int)*widthMS*heightMS); 
Rband->RasterIO(GF_Read,0,0,widthMS,heightMS,Rbuffer,widthMS,heightMS,GDT_Int32,0,0); 
Gband->RasterIO(GF_Read,0,0,widthMS,heightMS,Gbuffer,widthMS,heightMS,GDT_Int32,0,0); 
Bband->RasterIO(GF_Read,0,0,widthMS,heightMS,Bbuffer,widthMS,heightMS,GDT_Int32,0,0); 
Alphaband->RasterIO(GF_Read,0,0,widthMS,heightMS,Alphabuffer,widthMS,heightMS,GDT_Int32,0,0); 
MSimage = gcnew System::Drawing::Bitmap(widthMS,heightMS,PixelFormat::Format32bppArgb); 
k=0; 
for (j=0;j<heightMS;j++)
{
	for(i=0;i<widthMS;i++)
	{
		Color newclr = Color.FromArgb((unsigned char)(*(Alphabuffer+k)),(unsigned char)(*(Rbuffer+k)),(unsigned char)(*(Gbuffer+k)),(unsigned char)(*(Bbuffer+k)));
		MSimage->SetPixel(i,j,newclr); 
		k++; 
	}
} 
pictureBox2->Image = MSimage;

Posted 1-Dec-12 5:27am

Richa Vats

Updated 2-Dec-12 18:00pm

v3

Add a Solution

Comments

Sergey Alexandrovich Kryukov 1-Dec-12 22:19pm

I don't think anyone could help you unless you provide some code sample.
--SA

Richa Vats 2-Dec-12 8:12am

below is the code which i used to read a quickbird multispectral image.....

GDALDataset *hd1;
GDALRasterBand *Rband, *Gband, *Bband, *Alphaband;
int widthMS, heightMS;
int *Rbuffer, *Gbuffer, *Bbuffer, *Alphabuffer;
int i,j,k;
hd1 =(GDALDataset *)GDALOpen("G:\\images\\QB-multispectral-clipd.tif",GA_ReadOnly);
widthMS = hd1->GetRasterXSize();
heightMS = hd1->GetRasterYSize();
Rband = hd1->GetRasterBand(1);
Gband = hd1->GetRasterBand(2);
Bband = hd1->GetRasterBand(3);
Alphaband = hd1->GetRasterBand(4);
Rbuffer = (int*)CPLMalloc(sizeof(int)*widthMS*heightMS);
Gbuffer = (int*)CPLMalloc(sizeof(int)*widthMS*heightMS);
Bbuffer = (int*)CPLMalloc(sizeof(int)*widthMS*heightMS);
Alphabuffer = (int*)CPLMalloc(sizeof(int)*widthMS*heightMS);
Rband->RasterIO(GF_Read,0,0,widthMS,heightMS,Rbuffer,widthMS,heightMS,GDT_Int32,0,0);
Gband->RasterIO(GF_Read,0,0,widthMS,heightMS,Gbuffer,widthMS,heightMS,GDT_Int32,0,0);
Bband->RasterIO(GF_Read,0,0,widthMS,heightMS,Bbuffer,widthMS,heightMS,GDT_Int32,0,0);
Alphaband->RasterIO(GF_Read,0,0,widthMS,heightMS,Alphabuffer,widthMS,heightMS,GDT_Int32,0,0);

MSimage = gcnew System::Drawing::Bitmap(widthMS,heightMS,PixelFormat::Format48bppRgb);
k=0;
for (j=0;j<heightms;j++)
{
for="" (i="0;i<widthMS;i++)
" {
="" color="" newclr="Color::FromArgb(*(Rbuffer+k),*(Gbuffer+k),*(Bbuffer+k),*" (alphabuffer+k));
="" msimage-="">SetPixel(i,j,newclr);
k++;
}
}
pictureBox2->Image = MSimage;

Sergey Alexandrovich Kryukov 2-Dec-12 13:47pm

To make it readable, could you move it in the body of the question, using "Improve question" above?
To format code property, you need to inclose it in the HTML tag "pre" with the attribute lang="cs". But this does not work in comments.
--SA

Add your solution here

Treat my content as plain text, not as HTML

Preview 0

…

Existing Members

Sign in to your account

...or Join us

Download, Vote, Comment, Publish.

Your Email
Password
Forgot your password?

Your Email
This email is in use. Do you need your password?
Optional Password

I have read and agree to the Terms of Service and Privacy Policy
Please subscribe me to the CodeProject newsletters

When answering a question please:

Read the question carefully.
Understand that English isn't everyone's first language so be lenient of bad spelling and grammar.
If a question is poorly phrased then either ask for clarification, ignore it, or edit the question and fix the problem. Insults are not welcome.
Don't tell someone to read the manual. Chances are they have and don't get it. Provide an answer or move on to the next question.

Let's work to help developers, not make them feel stupid.

This content, along with any associated source code and files, is licensed under The Code Project Open License (CPOL)