Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1394 1747 (N41W56 L=293)   613    112        1 (1364/1699-1436/1774)  12.1 W49-W63 N38-N42
      2  491 1391 (N40E32 L=204)   666    118        0 ( 437/1360- 543/1417) -24.4 E38-E27 N37-N42
      3  832  962 (N07E13 L=224)   277     65        1 ( 814/ 921- 855/ 996)  78.7 E13-E11 N04-N09
      4 1812  837 (S23W53 L=289)  1526    166        6 (1766/ 763-1857/ 893)  26.6 W47-W58 S27-S20
      5 1693  615 (S35W41 L=277)   555    129        0 (1650/ 560-1726/ 667)  35.5 W36-W43 S38-S32
      6  743  601 (S13E26 L=210)   504    160        3 ( 709/ 540- 790/ 647)  50.2 E29-E22 S16-S10
      7  605  605 (S10E35 L=201)   345     69        0 ( 564/ 583- 638/ 620) -21.8 E37-E32 S11-S08
      8  759  486 (S21E29 L=207)   258     89        2 ( 715/ 470- 800/ 506)   0.0 E32-E25 S21-S19
      9 1485  400 (S46W21 L=257)   351     74        0 (1458/ 364-1505/ 429)  51.3 W18-W22 S48-S43
  

    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