Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  895 1826 (N39E39 L= 47)   406     84        2 ( 877/1798- 915/1851) -68.2 E41-E36 N36-N41
      2 1189 1649 (N36E10 L= 76)   863    132        2 (1135/1586-1236/1696) -26.6 E14-E04 N32-N38
      3  275 1520 (N04E71 L= 15)  1359    270        6 ( 235/1469- 298/1582) -60.3 E77-E65 N01-N08
      4  849 1413 (N11E21 L= 65)   311     71        2 ( 815/1383- 894/1433)  36.9 E22-E18 N08-N12
      5 1566 1278 (N24W25 L=111)   247     80        1 (1556/1233-1577/1320) -59.0 W23-W26 N21-N26
      6 1793 1058 (N19W48 L=134)   640     99        2 (1777/ 998-1813/1108) -74.1 W45-W50 N15-N22
      7 1583  920 (N04W35 L=121)   271     96        0 (1543/ 900-1614/ 954)  45.0 W32-W36 N02-N07
      8  479  773 (S35E28 L= 58)   831    325        2 ( 360/ 746- 621/ 801)  15.9 E38-E17 S38-S30
      9 1044  742 (S21W09 L= 95)   669    281        2 ( 910/ 701-1185/ 793)  16.4 W01-W18 S24-S18
  

    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