Проблема с преобразованием HSI-> RGB

У меня возникла странная проблема при попытке применить преобразование из HSI в RGB, которое я пытаюсь выполнить как часть домашнего задания. Для начала вот мое оригинальное изображение:

исходное изображение (извините за ссылки, у меня нет представителя для публикации изображений)

Я конвертирую из изображения со значениями H, S, I (0-255). Вот отображение этого изображения в цветовом пространстве RGB (оттенок = красный, насыщенный = зеленый, интенсивный = синий):

HSI SRC изображение

Для меня это выглядит правильно … белые цветы (высокой интенсивности) голубые, оттенки (красные), кажется, меняются с цветами разных цветов.

Так что теперь я делаю некоторое выравнивание на канале интенсивности (синий на дисплее) в некоторых полях, что приводит к следующему:

HSI выровненное изображение:
http://i.imgur.com/hgk9K.png
((Вы должны скопировать их, потому что у меня нет представителя для> 2 гиперссылок))

Ничего страшного там нет. Это выглядит не совсем правильно, но трудно сказать, и все, что находится за пределами области интересов, одинаково.

Теперь здесь возникает проблема. После преобразования этого обратно в RGB я получаю такой результат:

преобразованное изображение:
http://i.imgur.com/6wEyw.png

Понятия не имею, в чем здесь проблема, зеленый выглядит хорошо, но некоторые красные / синие пиксели выглядят максимально. Три поля там, как и ожидалось, и их содержимое кажется довольно запутанным, но это не самая большая проблема, так как это может быть из моей функции выравнивания. Все, что находится вне этих полей, теперь должно быть почти идентично исходному изображению, но каким-то образом оно перемешано.

Я много раз просматривал свой код и добавил дополнительные скобки и приведение типов данных просто для уверенности, но все еще не могу найти проблему. Я уверен, что мои формулы верны, но я, должно быть, упускаю некоторую проблему в способе вычисления значений пикселей.

Вот код, который я использую для преобразования из HSI в RGB. Возможно, ошибка не в этом методе, но все функции, использованные здесь, были протестированы в другом месте и, кажется, работают нормально.

void colorSpace::HSItoRGB(image &src, image &tgt){
cout<<"HSI->RGB\n";
tgt.resize(src.getNumberOfRows(),src.getNumberOfColumns());
float pi = 3.14159265358979f;
for (int i=0; i<src.getNumberOfRows(); i++){    //for each pixel
for (int j=0; j<src.getNumberOfColumns(); j++){
//re-normalize h,s,i
float h = ((float)src.getPixel(i,j,H))*pi*2.0f/255.0f;//255/2 instead of 180
float s = ((float)src.getPixel(i,j,S))/255.0f;//255 instead of 100
float in= ((float)src.getPixel(i,j,I))/255.0f;
//compute x y z
float x = in*(1.0f-s);
float y = in*( 1.0f + (s*cos(h) / cos(pi/3.0f-h)) );
float z = 3.0f*in-(x+y);
float r,g,b;    //set rgb
if(h<(2.0f*pi/3.0f)){
b = x;
r = y;
g = z;
}else if(h<(4.0f*pi/3.0f)){//&&2pi/3<=h
r = x;
g = y;
b = z;
}else{  //less than 2pi && 4pi/3<=h
g = x;
b = y;
r = z;
}
//convert normalized rgb to 0-255 range
int rr = (int)round(r*255.0f);
int gg = (int)round(g*255.0f);
int bb = (int)round(b*255.0f);
tgt.setPixel(i,j,RED,rr);
tgt.setPixel(i,j,GREEN,gg);
tgt.setPixel(i,j,BLUE,bb);
}
}
}

Кто-нибудь видит какие-либо проблемы в коде или есть понимание взгляда на изображения?