Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  593 1526 (N30E38 L=212)  1628    374        1 ( 526/1411- 697/1630) -26.6 E48-E26 N22-N36
      2 1594 1487 (N36W42 L=291)  1320    256        3 (1518/1368-1663/1566) -61.7 W36-W46 N28-N41
      3 1819  930 (N01W59 L=309)   674     74        2 (1790/ 902-1844/ 967)  36.9 W56-W62 N00-N03
      4 1488  927 (S01W30 L=280)   665    151        1 (1424/ 906-1568/ 961)   5.7 W26-W36 S02-N01
      5 1748  853 (S04W53 L=302)   616     65        1 (1721/ 817-1774/ 882) -76.0 W50-W55 S05-S01
      6 1818  776 (S08W63 L=312)   547     60        1 (1797/ 744-1838/ 812)  63.4 W60-W65 S10-S05
      7 1417  723 (S14W28 L=277)   985    279        2 (1307/ 657-1519/ 777) -19.4 W20-W37 S17-S10
      8 1221  687 (S18W16 L=265)   493    130        1 (1163/ 658-1270/ 717)   8.1 W11-W18 S19-S15
      9  621  657 (S25E25 L=224)   247     77        0 ( 578/ 644- 658/ 667)   0.0 E27-E21 S25-S24
     10 1372  545 (S26W29 L=278)   257     64        0 (1356/ 517-1388/ 582) -76.0 W27-W30 S27-S23
     11 1439  525 (S27W35 L=284)   292     90        0 (1410/ 493-1461/ 567)  80.5 W33-W36 S29-S23
  

    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