Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1253 1558 (N30W19 L=168)   808    201        0 (1176/1496-1337/1614) -41.6 W14-W23 N25-N34
      2  643 1497 (N29E24 L=126)   982    236        2 ( 562/1459- 749/1573) -31.0 E30-E15 N25-N34
      3  413 1348 (N21E40 L=110)  2174    333       10 ( 364/1226- 494/1478) -70.6 E45-E32 N12-N29
      4 1660 1451 (N21W47 L=197)   454    104        1 (1619/1410-1707/1473)  16.7 W42-W52 N18-N22
      5  269 1255 (N16E52 L= 98)   512     81        2 ( 240/1220- 309/1284)  38.7 E54-E48 N13-N17
      6 1440 1129 (N02W26 L=175)   472    178        2 (1369/1106-1526/1155) -11.3 W21-W31 N00-N03
      7 1580 1074 (S02W35 L=185)   414    147        2 (1541/ 999-1602/1123) -80.5 W32-W36 S05-N01
      8 1447  997 (S06W26 L=175)   261    116        0 (1400/ 951-1478/1025) -38.7 W22-W27 S08-S03
      9  196  846 (S08E61 L= 89)  3055    380       11 ( 123/ 713- 274/ 978) -55.5 E69-E53 S16-N00
     10 1407  891 (S12W23 L=173)   528     96        1 (1370/ 858-1436/ 937)  63.4 W20-W24 S13-S09
     11 1178  824 (S15W08 L=158)   996    357        3 (1069/ 722-1305/ 899)  32.7 W01-W16 S20-S10
     12 1331  849 (S14W18 L=168)   307     86        1 (1313/ 814-1352/ 891)  78.7 W16-W19 S16-S11
     13 1814  783 (S19W57 L=207)   464     87        0 (1783/ 753-1845/ 815)  33.7 W54-W60 S21-S17
     14  214  697 (S17E64 L= 86)  1239    166        1 ( 175/ 652- 257/ 738) -45.0 E68-E59 S19-S14
     15  512  512 (S31E41 L=108)   519     99        2 ( 463/ 486- 559/ 542)   6.3 E46-E37 S32-S28
     16 1696  455 (S41W59 L=208)  1300    178        1 (1608/ 432-1774/ 481)   0.0 W46-W73 S42-S39
     17  416  399 (S38E58 L= 91)  1331    170        3 ( 397/ 300- 439/ 466)  26.6 E75-E50 S44-S33
  

    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