Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1671 1501 (N32W53 L=108)  1087    166        1 (1640/1463-1713/1544) -80.5 W49-W56 N29-N35
      2 1461 1508 (N34W34 L= 89)   434     52        0 (1431/1475-1485/1541) -78.7 W32-W36 N31-N36
      3  323 1485 (N32E57 L=358)  1146    291        4 ( 287/1440- 383/1537) -28.6 E62-E49 N29-N35
      4  582 1391 (N27E31 L= 24)   249     63        1 ( 551/1378- 616/1407)   0.0 E33-E28 N26-N28
      5  742 1218 (N16E18 L= 37)   465    191        1 ( 694/1150- 801/1299) -51.3 E21-E13 N11-N21
      6  177 1243 (N16E67 L=348)  1264    197        2 ( 150/1193- 203/1295) -45.0 E71-E62 N13-N19
      7 1782 1207 (N12W56 L=111)   450     77        2 (1756/1171-1806/1233) -63.4 W53-W58 N10-N14
      8 1168  971 (N00W09 L= 64)   296     84        1 (1123/ 945-1212/ 995)   9.5 W05-W11 S01-N01
      9 1854  841 (S11W63 L=118)   650    130        1 (1842/ 800-1872/ 895) -63.4 W61-W66 S13-S07
     10  573  749 (S13E30 L= 25)   362    107        0 ( 519/ 723- 615/ 779) -26.6 E33-E27 S14-S11
     11  378  606 (S23E49 L=  6)  1628    422        3 ( 293/ 505- 455/ 701)  38.3 E60-E40 S31-S16
     12 1057  570 (S25W02 L= 57)   211     47        0 (1026/ 554-1087/ 584)   0.0 W00-W03 S26-S24
  

    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