IDL方程的解决方案

如何在IDL中用newton方法得到方程的解，我的程序将提供各种（在这里的方块），但牛顿方法ID IDL只是收回方程的初始解。请帮助我的论文进行遥感图像处理。IDL方程的解决方案

我的IDL程序是：

Pro TSM_lixiaModel 
!Except=0 
Compile_opt idl2 
dir='I:\lwkDATA\waterRegion\MODIS\' 
files=file_search(dir,'*RC.tif',count=num) 
for i=0,num-1 do begin 
file=files[i] 
    raster=e.Openraster(file,external_type='ADS40') 
    outFile=file_dirname(file)+'\'+file_basename(file,'.tif')$ 
    +'_TSM_lixiaModel.tif' 


TSMraster=enviraster(uri=outfile,$ 
     nrows=raster.nrows,$ 
     ncolumns=raster.ncolumns,$ 
     nbands=1,$ 
     data_type=raster.data_type,$ 
     SPATIALREF=raster.SPATIALREF) 
tileIterator=raster.createtileiterator(bands=2,tile_size=[100,100]) 
count=0 
foreach tile,tileiterator do begin 
    count++ 

    ***;Tile is variable, but not the needed Solution of the equation 
    ;S is needed Solution of the equation 
    ;????? is the initially solution of the ‘newtfunc’，how do i give the  various（Tile） to newtfunc*** 

    processedTile=newton(??????,'newtfunc'); 

    currentSubRect=tileIterator.Current_subrect 
    TSMraster.setdata,processedTile,sub_rect=currentsubrect 
    print,'1' 
endforeach 
TSMraster.save 
endfor 
End 

function newtfunc,S 
    compile_opt idl2 
    return,(r-93.0943)*S+49.2464*S*exp(-0.0001*S)-344.016+45*r 
end

来源

2017-03-01 wenkai li

我不完全理解你的代码，但也许你只需要使用上ENVIRaster GetData方法？例如：

data = raster.GetData(BANDS=[0], SUB_RECT=[100,449,550,899])

然后，您可以将您的数据传递到牛顿函数。完整的文档是在这里： https://www.harrisgeospatial.com/docs/enviraster__getdata.html

顺便说一句，在一个相关的说明，如果你需要将其他参数传递到牛顿的功能，一招是使用公共块，在这里建立了共同的在你的调用例程中阻塞，然后访问你的“newtfunc”中的变量。例如：

pro mymain 
    common foo, x, y, z 
    x=1 & y=1 & z=1 
    result = newton(...,'newtfunc') 
end 

function newtfunc,S 
    common foo 
    <use x,y,z here> 
    return,... 
end

希望这可以帮助一些！如果不是，您可以联系技术支持。

来源

2017-03-01 22:44:18

谢谢。函数包括R和S，牛顿方法刚开始需要S，那么牛顿方法就会计算出真实的S.但是在这里，我需要先给出S和各种R，R需要用牛顿法来计算。那我的麻烦 –

IDL方程的解决方案

回答

相关问题