Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1485 1696 (N46W30 L=210)   404     65        0 (1460/1663-1510/1721)  36.9 W26-W34 N42-N48
      2  733 1585 (N26E27 L=154)   275     25        0 ( 707/1566- 756/1604)  45.0 E28-E24 N25-N27
      3 1780 1393 (N28W53 L=233)   380     42        0 (1759/1366-1799/1416) -33.7 W50-W55 N26-N29
      4  296  928 (S17E50 L=130)   437     25        0 ( 273/ 908- 325/ 950) -14.0 E51-E46 S18-S15
      5  909  692 (S27E04 L=176)   244     78        0 ( 873/ 659- 932/ 729)  59.0 E05-E02 S29-S24
      6 1098  697 (S24W08 L=189)   354    148        0 (1022/ 669-1163/ 717)   0.0 W03-W13 S25-S23
      7 1719  639 (S19W56 L=236)   524     58        1 (1692/ 605-1740/ 670) -26.6 W52-W59 S20-S16
      8  551  642 (S33E30 L=150)   256     38        0 ( 537/ 614- 569/ 666)  63.4 E31-E28 S35-S31
      9 1595  476 (S31W53 L=233)   599    142        0 (1556/ 418-1650/ 528)  82.9 W50-W55 S34-S26
     10 1013  259 (S57W14 L=194)   737     95        2 ( 965/ 228-1055/ 290)   0.0 W08-W18 S59-S54
     11 1123  239 (S57W27 L=207)   798     75        2 (1093/ 209-1162/ 278) -22.6 W21-W33 S59-S53
  

    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