Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  534 1713 (N45E51 L=191)   633    112        1 ( 502/1687- 573/1740)   0.0 E54-E46 N42-N46
      2 1598 1587 (N43W46 L=288)   632    150        1 (1554/1563-1648/1610)  10.3 W40-W51 N40-N44
      3 1328 1360 (N25W19 L=261)   244     83        1 (1297/1342-1369/1371)  21.8 W16-W21 N23-N25
      4 1617 1177 (N15W39 L=281)   387     84        1 (1587/1153-1645/1220) -53.1 W36-W41 N13-N17
      5 1798 1071 (N09W56 L=298)  1034    205        1 (1765/ 969-1822/1161) -76.0 W52-W58 N02-N14
      6 1877 1013 (N05W66 L=308)   767    143        2 (1858/ 967-1901/1068) -63.4 W63-W70 N02-N08
      7 1001  887 (S06E01 L=241)   793    214        2 ( 957/ 772-1025/ 986) -72.9 E03-W01 S13-N00
      8  551  671 (S23E31 L=211)   252     88        0 ( 511/ 657- 591/ 684)  18.4 E33-E27 S24-S21
      9 1416  597 (S23W30 L=271)   263     73        0 (1385/ 583-1451/ 619) -11.3 W27-W32 S23-S21
     10  365  531 (S35E50 L=192)  1819    451        2 ( 254/ 508- 448/ 568)  12.3 E64-E41 S38-S31
  

    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