Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1620 1666 (N32W62 L= 30)   585    180        0 (1592/1636-1636/1701)  14.0 W57-W65 N31-N34
      2 1409 1540 (N30W38 L=  6)  1153    246        3 (1320/1497-1502/1587)  -3.8 W31-W46 N28-N32
      3  486 1483 (N44E35 L=293)   377     74        1 ( 448/1468- 522/1497)   0.0 E38-E31 N42-N44
      4  808 1206 (N21E11 L=317)   215     95        0 ( 754/1179- 851/1230) -38.7 E13-E08 N18-N22
      5 1555  937 (S08W32 L=360)  1303    272        3 (1506/ 819-1589/1070)  63.4 W27-W35 S14-N00
      6 1777 1017 (S08W51 L= 19)   690    170        2 (1733/ 958-1818/1067) -53.1 W48-W54 S12-S04
      7 1683  939 (S11W42 L= 10)   718    181        2 (1657/ 854-1708/1034)  59.0 W38-W45 S15-S05
      8  222  920 (N10E60 L=268)  1053    215        2 ( 172/ 869- 279/ 987)  23.2 E67-E52 N06-N13
      9 1845  894 (S18W57 L= 25)   557    141        2 (1832/ 849-1859/ 940)  59.0 W55-W58 S20-S15
     10 1303  434 (S37W09 L=337)  1284    294        0 (1211/ 351-1366/ 559)  42.9 W00-W14 S42-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