Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  725 1636 (N20E38 L=281)   486    124        1 ( 698/1568- 751/1675) -33.7 E41-E32 N16-N22
      2 1521 1318 (N24W21 L=340)   311    116        0 (1510/1266-1533/1380) -66.8 W19-W23 N21-N28
      3 1769 1232 (N28W43 L=  2)  1708    348        2 (1726/1102-1819/1360) -68.2 W38-W49 N19-N35
      4 1856 1097 (N24W55 L= 14)   395     92        0 (1848/1050-1865/1143) -80.5 W53-W56 N20-N26
      5  139 1015 (S28E70 L=248)  1266    106        2 (  80/ 991- 183/1038) -90.0 E89-W00 S29-S26
      6  674  927 (S21E18 L=301)   247     85        0 ( 653/ 881- 702/ 959) -78.7 E18-E16 S24-S19
      7  958  731 (S25W05 L=323)   330    162        0 ( 889/ 711-1039/ 755)  40.6 W01-W08 S27-S21
      8 1715  649 (S05W56 L= 15)   325     49        2 (1696/ 622-1736/ 669)  90.0 W54-W58 S07-S03
      9 1404  483 (S24W42 L=  1)   974    233        1 (1299/ 447-1516/ 518)  42.3 W37-W48 S29-S19
     10 1004  393 (S44W23 L=341)   256     54        0 ( 969/ 378-1035/ 408)   9.5 W19-W25 S44-S42
     11  551  388 (S62E18 L=300)   420     83        0 ( 505/ 375- 591/ 401)  13.0 E25-E12 S64-S59
  

    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