Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1110 1890 (N57E35 L=  1)   885     67        1 (1084/1861-1136/1919) -17.4 E44-E28 N53-N60
      2 1476 1716 (N56W12 L= 49)   647     83        1 (1449/1667-1504/1745)  77.5 W09-W16 N51-N61
      3 1683 1580 (N53W40 L= 76)   690     76        3 (1660/1539-1706/1615)  81.9 W36-W43 N49-N56
      4 1616 1258 (N26W31 L= 67)   355    124        0 (1564/1225-1642/1307)  60.3 W27-W32 N22-N29
      5 1627 1106 (N17W34 L= 71)   595    171        1 (1601/1022-1670/1168)  90.0 W33-W37 N11-N21
      6 1761 1050 (N18W47 L= 84)   482     57        0 (1732/1024-1790/1072)  -9.5 W44-W50 N16-N19
      7 1892  974 (N18W65 L=101)  1916    145        4 (1853/ 922-1933/1031)  58.0 W59-W74 N14-N22
      8  238  934 (S29E55 L=341)   961    205        2 ( 201/ 821- 268/1007) -70.0 E60-E51 S35-S24
      9 1708  965 (N12W44 L= 81)   435     81        1 (1659/ 936-1748/ 983)  20.6 W39-W47 N10-N13
     10 1824  902 (N12W57 L= 94)  1017    123        1 (1786/ 854-1852/ 972)  76.0 W54-W61 N08-N16
     11  356  865 (S31E42 L=355)   436    131        0 ( 299/ 835- 414/ 891)  24.4 E47-E36 S33-S27
     12  186  834 (S37E64 L=333)   687     74        3 ( 157/ 798- 219/ 867) -76.0 E69-E58 S38-S34
     13  447  787 (S33E32 L=  4)   372    115        0 ( 396/ 769- 497/ 823)  29.1 E37-E28 S35-S29
     14 1191  414 (S35W32 L= 69)   623     80        2 (1151/ 395-1239/ 443)  31.0 W28-W35 S37-S32
     15  979  419 (S41W18 L= 54)   308    101        1 ( 923/ 404-1030/ 438)  20.6 W13-W21 S43-S39
     16  502  392 (S62E27 L=  9)   390     55        0 ( 476/ 369- 529/ 410)  -8.1 E31-E23 S63-S60
  

    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