Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  717 1742 (N62W04 L= 19)   505     59        1 ( 693/1713- 743/1767)  14.0 W00-W08 N58-N64
      2 1127 1702 (N42W33 L= 48)   739    153        1 (1064/1659-1194/1742) -77.5 W29-W35 N37-N46
      3  667 1690 (N59E04 L= 11)   333     59        1 ( 638/1666- 689/1718) -38.7 E07-E01 N56-N61
      4 1483 1682 (N28W57 L= 71)   674    153        0 (1398/1668-1554/1694) -24.8 W49-W62 N25-N31
      5  568 1632 (N58E17 L=357)   706    111        3 ( 522/1591- 616/1668)   0.0 E23-E10 N55-N60
      6  351 1520 (N56E47 L=327)   425     49        1 ( 333/1497- 369/1549)  56.3 E50-E43 N53-N57
      7 1171 1518 (N30W25 L= 39)   302     19        2 (1147/1497-1194/1540) -18.4 W23-W26 N27-N30
      8 1234 1337 (N17W21 L= 36)   426    221        1 (1133/1316-1341/1370) -36.0 W15-W26 N13-N22
      9  302 1118 (N29E48 L=327)   307     66        0 ( 268/1099- 331/1142)   8.1 E51-E44 N28-N30
     10  419  980 (N18E39 L=336)   234     76        1 ( 391/ 953- 439/1015)  21.8 E41-E36 N17-N20
     11 1465  622 (S30W16 L= 30)   635    169        0 (1397/ 585-1514/ 681)  11.3 W10-W20 S31-S27
     12 1247  322 (S44E09 L=  6)   999    187        2 (1161/ 272-1320/ 374)   7.1 E17-E01 S46-S41
     13  898  222 (S40E40 L=335)   723    131        1 ( 853/ 176- 934/ 278)  65.6 E43-E36 S45-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