Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  192 1159 (S11E62 L=262)  1086    240        1 ( 130/1086- 244/1243)  48.0 E68-E56 S15-S05
      2 1770 1049 (N11W48 L= 12)   717    190        1 (1726/ 959-1794/1123) -53.1 W43-W50 N05-N14
      3 1863  995 (N10W58 L= 22)   846    268        3 (1826/ 900-1904/1121) -50.2 W52-W66 N05-N17
      4  333  985 (S20E45 L=278)  2188    683        0 ( 154/ 836- 558/1078)  25.9 E62-E28 S30-S14
      5 1904  866 (N04W68 L= 32)  1449    218        0 (1868/ 824-1943/ 898)  33.7 W64-W74 N01-N07
      6  829  739 (S27E07 L=317)   278     54        0 ( 802/ 710- 858/ 774)  68.2 E08-E05 S29-S24
      7  361  643 (S41E43 L=280)  1074    233        1 ( 306/ 573- 414/ 714)  42.5 E50-E37 S46-S35
      8 1493  648 (S18W38 L=  1)   461    156        0 (1417/ 630-1549/ 690)  36.9 W32-W40 S20-S14
      9 1751  649 (S12W59 L= 22)   578    177        0 (1695/ 619-1793/ 667)   4.8 W52-W64 S13-S10
     10 1809  593 (S13W68 L= 32)   634    108        1 (1778/ 571-1836/ 623)  38.7 W64-W71 S14-S10
     11 1558  402 (S31W54 L= 18)  3082    416        5 (1421/ 367-1738/ 447)  28.4 W43-W69 S35-S22
  

    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