Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1011 1764 (N38E26 L=100)   944    181        2 ( 956/1707-1100/1791)  68.2 E28-E21 N31-N42
      2  448 1622 (N14E57 L= 69)   390     46        2 ( 423/1601- 476/1640)  36.9 E59-E54 N12-N15
      3 1468 1500 (N34W14 L=140)   370     26        0 (1449/1473-1487/1527) -33.7 W12-W15 N31-N35
      4  919 1344 (N08E14 L=112)   585    190        2 ( 827/1325-1011/1368)  26.6 E19-E09 N05-N10
      5 1655 1308 (N28W32 L=158)   461    120        0 (1623/1258-1686/1357) -29.7 W28-W35 N25-N29
      6  451 1218 (S10E40 L= 87)   309    136        0 ( 412/1160- 482/1279)  82.9 E41-E38 S13-S05
      7  870 1111 (S06E11 L=115)   441    147        0 ( 825/1058- 907/1176) -35.5 E14-E07 S08-S03
      8  610 1119 (S12E26 L=100)   302     70        0 ( 590/1090- 639/1154)  76.0 E27-E24 S13-S09
      9  847 1039 (S11E10 L=116)   443    186        1 ( 804/ 980- 898/1113) -29.1 E14-E05 S14-S07
     10 1893 1037 (N20W60 L=186)  1307    180        3 (1848/ 984-1940/1104) -20.0 W53-W69 N17-N23
     11 1799  961 (N12W51 L=177)   606    164        1 (1746/ 923-1848/ 989)  33.7 W47-W56 N08-N14
     12  960  942 (S13E01 L=125)   210     91        1 ( 924/ 926-1014/ 962)  26.6 E03-W01 S14-S11
     13 1804  845 (N06W55 L=181)  1247    188        5 (1742/ 771-1861/ 915) -12.1 W47-W61 N03-N09
     14 1113  811 (S16W11 L=137)   237    110        2 (1067/ 784-1159/ 831)  31.0 W08-W13 S18-S14
     15 1073  569 (S31W17 L=143)   323     77        1 (1039/ 540-1106/ 594)  -9.5 W13-W19 S32-S29
     16 1667  499 (S16W61 L=187)   787     96        2 (1622/ 477-1712/ 521)  18.4 W56-W66 S17-S13
     17 1406  393 (S30W48 L=174)   488     80        3 (1365/ 370-1434/ 425)  90.0 W45-W51 S32-S27
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y