Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1290 1520 (N32W02 L=123)   316     73        0 (1260/1483-1321/1548) -36.9 W00-W05 N30-N33
      2 1709 1407 (N40W39 L=160)   420     68        0 (1683/1388-1741/1428)   0.0 W36-W42 N37-N40
      3 1632 1276 (N27W32 L=153)  1297    242        1 (1559/1172-1704/1333) -19.7 W25-W39 N22-N30
      4  552 1294 (S01E35 L= 86)   741    134        1 ( 484/1243- 595/1331) -11.3 E41-E31 S03-N00
      5 1500 1270 (N23W22 L=143)   443    117        1 (1466/1215-1523/1324)  82.9 W20-W23 N18-N26
      6 1391 1164 (N13W17 L=138)   303     71        1 (1350/1142-1430/1185)  45.0 W15-W19 N11-N15
      7  324 1135 (S15E49 L= 72)   374    139        0 ( 282/1097- 393/1158)  33.7 E52-E43 S18-S12
      8 1793 1018 (N17W51 L=172)  1050    150        2 (1742/ 959-1829/1084) -26.6 W44-W55 N14-N20
      9 1846 1047 (N21W57 L=178)   382     39        1 (1828/1027-1867/1065)  36.9 W54-W59 N19-N22
     10  172  943 (S30E65 L= 56)  2172    238        2 ( 133/ 865- 214/1054)  40.9 E74-E59 S35-S22
     11 1653  794 (N01W45 L=166)   857    211        2 (1624/ 713-1681/ 903) -74.7 W43-W47 S04-N07
     12 1344  795 (S08W25 L=145)  1677    427        5 (1152/ 735-1494/ 857)  37.9 W14-W32 S15-S01
     13  184  802 (S39E66 L= 55)   645     52        2 ( 156/ 779- 209/ 826)  18.4 E72-E60 S40-S36
     14  704  696 (S32E11 L=110)  1012    268        2 ( 621/ 600- 814/ 776)  56.3 E15-E04 S40-S24
     15 1704  553 (S10W62 L=183)   526     52        0 (1673/ 536-1734/ 571)  45.0 W59-W65 S12-S08
     16 1436  313 (S32W59 L=180)  2111    309        2 (1293/ 271-1561/ 342)  14.5 W45-W79 S37-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