Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  836 1822 (N62W20 L=118)   489    110        0 ( 800/1784- 881/1864)  15.5 W12-W30 N59-N64
      2 1141 1785 (N46W39 L=137)   528     88        1 (1112/1746-1176/1810)   0.0 W35-W43 N43-N48
      3 1399 1696 (N32W50 L=148)  1603    251        3 (1305/1603-1484/1754) -81.3 W44-W54 N24-N37
      4  932 1687 (N47W16 L=114)   284     49        1 ( 901/1666- 952/1710)   0.0 W13-W18 N45-N48
      5  771 1589 (N45E00 L= 98)   277     78        0 ( 744/1552- 792/1615) -80.5 E02-W00 N41-N47
      6  947 1350 (N25W04 L=102)   347     98        2 ( 920/1312- 976/1402)  80.5 W01-W05 N21-N27
      7  523 1384 (N38E25 L= 73)   314     94        0 ( 483/1359- 561/1401) -29.7 E28-E21 N35-N39
      8  244 1378 (N45E58 L= 40)   736     90        2 ( 213/1341- 271/1399) -18.4 E64-E52 N42-N46
      9 1103 1248 (N15W10 L=108)   270    130        2 (1067/1220-1147/1279) -11.3 W07-W12 N13-N16
     10  204 1080 (N26E58 L= 40)  1825    276        2 ( 165/ 998- 274/1197) -58.4 E64-E50 N20-N33
     11 1244 1081 (N02W14 L=112)   248     69        1 (1219/1054-1273/1118)  21.8 W11-W16 N00-N03
     12  253  923 (N16E55 L= 43)   442     77        1 ( 228/ 886- 282/ 960)  53.1 E58-E51 N14-N18
     13 1118  811 (S10E00 L= 98)   638    211        0 (1085/ 703-1155/ 895)  55.0 E04-W03 S15-S05
     14 1753  646 (S39W44 L=142)   873    122        2 (1700/ 618-1821/ 678)   0.0 W37-W53 S40-S37
     15  898  463 (S25E23 L= 75)   389    128        0 ( 833/ 447- 959/ 478) -20.6 E27-E19 S26-S23
     16 1243  240 (S53E12 L= 86)   684     67        1 (1215/ 211-1284/ 268)   7.1 E15-E07 S56-S50
     17 1122  166 (S56E29 L= 69)   458     50        0 (1097/ 148-1144/ 194)  12.5 E33-E24 S58-S53
  

    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