Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  926 1895 (N60E46 L=121)   634     55        0 ( 903/1877- 956/1919)  -7.1 E56-E39 N58-N62
      2 1499 1539 (N42W27 L=194)   604    112        3 (1461/1480-1531/1583) -76.0 W23-W29 N37-N45
      3  439 1499 (N18E51 L=116)   685     86        1 ( 407/1458- 486/1535) -16.7 E56-E45 N15-N20
      4  376 1424 (N12E54 L=113)   552     72        2 ( 350/1398- 413/1453)  18.4 E56-E50 N10-N14
      5 1795 1278 (N29W55 L=222)  1858    241        2 (1707/1144-1866/1422) -48.8 W46-W63 N21-N37
      6  583 1352 (N11E34 L=133)   245     91        0 ( 546/1329- 625/1375)  38.7 E36-E31 N09-N13
      7  190 1106 (S10E64 L=103)   866    110        2 ( 153/1065- 236/1140)  24.4 E69-E58 S12-S07
      8  275 1086 (S10E53 L=114)   689    100        1 ( 236/1043- 313/1121) -45.0 E57-E49 S12-S07
      9  150  905 (S23E70 L= 97)   795     85        0 ( 132/ 864- 176/ 936)  26.6 E75-E64 S25-S20
     10  253  841 (S25E55 L=112)   418     86        0 ( 222/ 821- 288/ 871) -14.0 E58-E50 S26-S23
     11  877  668 (S25E03 L=164)  2446    575        5 ( 640/ 559-1093/ 806)  42.6 E18-W08 S36-S12
     12 1357  639 (S18W29 L=196)   407    140        0 (1321/ 615-1395/ 684)  90.0 W26-W31 S19-S14
     13 1700  645 (S11W56 L=223)   493     63        3 (1669/ 620-1729/ 672)  45.0 W53-W58 S12-S09
     14 1593  435 (S26W58 L=225)   827    125        1 (1541/ 402-1660/ 468)  41.2 W53-W63 S29-S22
     15 1018  184 (S61W30 L=197)   680     67        2 ( 982/ 165-1058/ 204)  33.7 W24-W35 S63-S58
     16 1235  177 (S55W54 L=221)   548     51        1 (1215/ 157-1258/ 203) -12.1 W48-W62 S56-S51
  

    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