Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1365 1749 (N48W47 L=201)   393     74        0 (1336/1718-1395/1776)  14.0 W41-W53 N46-N49
      2 1492 1751 (N45W62 L=216)   514     56        0 (1460/1731-1518/1771)   4.4 W55-W68 N44-N46
      3  643 1694 (N56E29 L=125)   425     42        0 ( 621/1665- 663/1722) -78.7 E31-E27 N52-N58
      4  565 1354 (N30E29 L=125)  1313    280        5 ( 431/1326- 684/1374)   0.0 E40-E19 N28-N32
      5 1684 1232 (N09W48 L=202)   278     62        0 (1657/1205-1705/1256) -63.4 W45-W49 N07-N11
      6  488 1104 (N15E34 L=120)   539    190        0 ( 448/1033- 523/1184)  71.6 E37-E31 N11-N20
      7  921  972 (N04E07 L=147)  1050    389        1 ( 837/ 803- 966/1114) -78.1 E10-E04 S06-N13
      8 1614  957 (S06W37 L=192)   571     96        3 (1582/ 916-1648/1003) -56.3 W35-W39 S08-S02
      9 1454  893 (S07W25 L=180)   207    122        0 (1439/ 840-1467/ 952)  81.9 W24-W26 S10-S03
     10  768  761 (S08E19 L=136)   617    287        0 ( 660/ 728- 900/ 781)  -7.1 E26-E10 S10-S06
     11 1815  717 (S25W58 L=212)   956     95        3 (1793/ 670-1832/ 767) -74.1 W55-W60 S28-S21
     12 1748  668 (S27W51 L=205)   744    144        1 (1723/ 609-1775/ 734)  82.9 W48-W53 S30-S22
     13 1270  615 (S22W11 L=166)   928    159        3 (1214/ 551-1325/ 686)  53.1 W07-W15 S26-S18
     14 1431  485 (S34W23 L=177)  3082    512        4 (1267/ 396-1575/ 646)  31.0 W10-W34 S40-S24
     15  808  467 (S28E22 L=133)   459    132        0 ( 758/ 445- 862/ 494) -23.2 E24-E17 S29-S25
     16  609  396 (S31E40 L=114)   598    109        1 ( 572/ 361- 653/ 437)  24.0 E43-E34 S33-S28
     17 1448  278 (S54W32 L=186)  1066    136        2 (1373/ 256-1504/ 302) -10.0 W21-W40 S57-S50
  

    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