Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  702 1660 (N32E39 L=197)   603    130        0 ( 648/1621- 737/1685)   0.0 E43-E34 N29-N33
      2 1378 1646 (N47W13 L=249)   370     85        1 (1333/1622-1418/1662)   0.0 W08-W17 N45-N48
      3 1017 1569 (N32E13 L=223)   759    125        2 ( 969/1523-1069/1608) -20.6 E16-E08 N29-N34
      4  664 1496 (N20E34 L=202)   345     75        0 ( 624/1483- 707/1513)  21.8 E36-E31 N18-N21
      5 1577 1132 (N16W33 L=270)   557     90        2 (1554/1092-1603/1172)  78.7 W31-W35 N13-N18
      6  166 1065 (S15E67 L=169)  1167    124        2 ( 139/1012- 198/1111)  41.2 E71-E62 S18-S11
      7 1681 1064 (N14W43 L=279)   356     85        2 (1656/1013-1711/1098) -53.1 W40-W45 N11-N15
      8 1426  727 (S11W31 L=267)   755    162        1 (1352/ 698-1495/ 754)   5.7 W26-W36 S13-S09
      9 1248  661 (S18W22 L=258)   241     97        1 (1217/ 624-1284/ 692)  71.6 W20-W23 S21-S15
     10 1470  611 (S17W38 L=274)   769    258        3 (1358/ 574-1571/ 647)  31.6 W31-W45 S21-S13
     11 1491  285 (S37W65 L=302)   685     98        1 (1459/ 265-1534/ 311)  20.6 W58-W71 S38-S35
  

    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