Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1149 1697 (N37E14 L=298)   490     91        1 (1100/1670-1200/1718)  15.9 E18-E10 N35-N38
      2  444 1543 (N09E53 L=260)   549     77        0 ( 415/1509- 478/1577) -20.6 E56-E48 N07-N10
      3 1626 1514 (N42W27 L=339)   723    140        1 (1596/1452-1651/1568) -60.3 W23-W29 N38-N46
      4 1586 1285 (N24W26 L=339)   410    158        0 (1519/1236-1615/1321)   0.0 W21-W29 N21-N26
      5 1800 1133 (N23W47 L=360)   627    117        0 (1758/1072-1836/1171) -20.6 W42-W52 N20-N25
      6 1735  974 (N12W45 L=357)   379     19        1 (1709/ 956-1761/ 999)  21.8 W42-W47 N10-N13
      7  773  914 (S19E11 L=301)   229     60        0 ( 740/ 893- 807/ 932)  11.3 E13-E08 S19-S17
      8 1649  881 (N04W41 L=353)   243     90        0 (1598/ 870-1689/ 894)  26.6 W37-W43 N02-N05
      9 1278  715 (S16W23 L=336)   306     66        1 (1243/ 702-1326/ 732)  26.6 W21-W26 S17-S14
     10  471  581 (S48E28 L=285)   376     47        3 ( 445/ 553- 497/ 609) -11.3 E30-E24 S49-S45
     11 1509  388 (S25W56 L=  9)   530     94        0 (1471/ 358-1549/ 412)   0.0 W51-W62 S26-S24
     12  934  250 (S56W29 L=342)   418     62        1 ( 917/ 214- 950/ 287) -26.6 W24-W35 S58-S52
  

    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