Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1633 1543 (N15W57 L=128)   494    102        0 (1596/1497-1659/1586) -90.0 W55-W58 N12-N18
      2 1162 1370 (N21W18 L= 89)  1784    463        1 ( 968/1244-1327/1445) -56.3 W08-W23 N09-N30
      3 1696 1164 (S06W44 L=115)  1198    261        2 (1570/1105-1797/1205)  -8.1 W33-W54 S08-S03
      4  483 1051 (N21E32 L= 39)   240     67        0 ( 468/1025- 509/1083) -56.3 E33-E30 N19-N23
      5 1698  933 (S18W38 L=109)  1636    291        2 (1642/ 805-1743/1063)  47.5 W31-W45 S24-S12
      6 1781  896 (S24W45 L=116)   489     53        0 (1764/ 861-1797/ 932)  63.4 W43-W46 S25-S21
      7  974  630 (S14E14 L= 57)   565    184        2 ( 938/ 561- 995/ 723)  45.0 E17-E09 S17-S09
      8  849  655 (S10E20 L= 51)   370     74        2 ( 823/ 624- 873/ 694)  36.9 E21-E17 S11-S07
      9  750  480 (S18E32 L= 39)  1213    309        1 ( 642/ 372- 809/ 602) -90.0 E34-E30 S25-S07
     10  262  538 (S04E70 L=  0)  1177    103        1 ( 242/ 485- 285/ 590)  20.6 E80-E63 S07-S01
     11  508  469 (S13E49 L= 22)   940    175        2 ( 434/ 440- 580/ 493) -29.1 E56-E44 S15-S09
     12  844  188 (S41E47 L= 24)  1541    153        2 ( 810/ 113- 876/ 247)  34.4 E61-E40 S47-S34
  

    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