Problem with "Index out of range" error.

Share icon

Hello,

I am trying to create a code that firstly creates a dataset that simulates a radar measurements of a cloud in the atmosphere. Then I want to convert it into data with real height and range. However, this is not the problem. I get an "Index out of range" error for the RealCloud wave. When I debug me code it seems that it stucks in the RealCloudConv function, and it doesn't return the value into RCloud in fake2real function. Could anyone help me please?

This is my whole code (there are 4 different functions. The first is the basic. The two last are secondary functions and there is no problem with them. The problem in the last command of the first function or somewhere in the second function.

Function Fake2Real(HoRes,VeRes) // Calculates real height and real range of radar values in two different matrices

variable HoRes,VeRes

variable Perimeter, RadiusOfRange, RangeColumns, i, j, HeightRaws, Rad, EarthRadius, a, b, c, d

variable  RealHeightMax, RealRangeMax

make/O/N=(60/VeRes,200/HoRes) RealRange

make/O/N=(60/VeRes,200/HoRes) RealHeight

make/O/N=(60/VeRes,200/HoRes) cloud=0

wave cloud20

LoadWave/J/D/O/M/W/K=0/N=cloud2 "\\\\nask.man.ac.uk\\home$:Desktop:FakeCloud.csv"

RangeColumns = 200/HoRes  //200 is the total range of the radar

HeightRaws = 60/VeRes

Perimeter=44050

EarthRadius = 6371

Rad=Deg2Rad(VeRes)

    For (i=0; i<RangeColumns; i+=1)
       
        For (j=0; j<HeightRaws; j+=1)

        RadiusOfRange = ((i+1)*HoRes*360)/Perimeter

        RealRange[j][i]= (RadiusOfRange*2*pi*6371)/360
       
        RealHeight[j][i]= ((RealRange[j][i])*sin(Rad))+((RealRange[j][i]^2)/(2*EarthRadius))*(HeightRaws-j)
           
        EndFor
    EndFor
   
   
RealRangeMax=round(Maximum(HeightRaws,RangeColumns,RealRange))

RealHeightMax=round(Maximum(HeightRaws,RangeColumns,RealHeight))

Make/O/N=(RealHeightMax, RealRangeMax) RCloud=0

RCloud=RealCloudConv(cloud20, RealRange, RealHeight, VeRes, HoRes, RealHeightMax, RealRangeMax)


End

//-------------------------------------------------------------------------------------------------------------------------//3

Function RealCloudConv(cloud1, RealRange1, RealHeight1, VeRes1, HoRes1, RealHeightMax1, RealRangeMax1)

wave cloud1, RealRange1, RealHeight1

variable VeRes1, HoRes1, RealHeightMax1, RealRangeMax1

variable i, j, raws, columns, col, raw, realcols, RealRangeMax, RealHeightMax

raws=dimsize(cloud1,0)
columns=dimsize(cloud1,1)

make/O/N=(RealHeightMax1, RealRangeMax1) RealCloud=0 //60/VeRes1,200/HoRes1

    For (i=0; i<raws; i+=1)
       
        For (j=0; j<columns; j+=1)
       
            col=round(RealRange1[i][j])
            raw=round(RealHeight1[i][j])
       
            RealCloud[raw][col] = cloud1[i][j]     
       
        EndFor 
   
    EndFor
   
return RealCloud
   
End


//--------------------------------------------------------------------------------------------------------------------------//
//SECONDARY FUNCTIONS//
//--------------------------------------------------------------------------------------------------------------------------//

Function Maximum(raws1,columns1,matrix)

wave matrix

variable raws1,columns1

variable max1, i, j

max1=0

for (i=0; i<raws1; i+=1)
    for(j=0; j<columns1; j+=1)
        if (matrix[i][j]>max1)
            max1=matrix[i][j]
        endif
    endfor
endfor

return max1

end

//------------------------------------------------------------------------------------------------------------------------//

Function Deg2Rad(VRes)   //Converts degrees to radians

variable VRes

variable radians

radians=VRes*0.005555556*pi

return radians

End
Could you maybe use the igor tags and cut the line length to 120 character?
This will make the code much easier to read.

One mistake is that RealCloudConv returns a wave but declares the return type to be a variable.

May 5, 2015 at 03:43 am - Permalink

Hello and thank you for your reply. I got it. You mean Function/wave.

Thanks a lot.

May 5, 2015 at 07:40 am - Permalink

Functions return a numeric variable unless otherwise specified. Have a look at
displayhelptopic "Function"

Igor tags are "igor" and "/igor" each enclosed in "<" and ">".
HJ

May 5, 2015 at 07:44 am - Permalink