Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1624 1543 (N20W58 L=187)   632    114        0 (1573/1524-1684/1562) -32.0 W53-W62 N17-N22
      2 1250 1489 (N29W27 L=156)   389    101        1 (1217/1458-1279/1529)  26.6 W24-W30 N26-N30
      3 1628 1464 (N16W53 L=182)   383     45        0 (1604/1437-1650/1485) -71.6 W51-W54 N14-N17
      4  193 1047 (N23E63 L= 66)  1824    215        5 ( 146/ 954- 237/1117) -45.0 E71-E57 N17-N27
      5 1847  942 (S20W57 L=186)   539     55        2 (1828/ 915-1872/ 972)  33.7 W54-W60 S21-S17
      6 1624  858 (S17W34 L=163)   687    164        2 (1578/ 794-1680/ 918)  29.1 W29-W39 S20-S14
      7  826  489 (S21E25 L=104)  1748    416        3 ( 698/ 356- 931/ 638)  27.5 E39-E14 S27-S12
      8  594  367 (S25E48 L= 81)  1849    378        2 ( 527/ 296- 701/ 508)  79.4 E54-E41 S31-S14
      9  938  460 (S25E20 L=110)   561    162        1 ( 880/ 414- 997/ 494) -66.8 E23-E16 S29-S22
     10 1075  437 (S30E12 L=118)   331     77        0 (1032/ 416-1117/ 465)   0.0 E14-E08 S30-S28
  

    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