Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  642 1833 (N34E61 L=264)   843     80        1 ( 610/1804- 676/1870) -71.6 E65-E57 N30-N37
      2  989 1741 (N35E26 L=299)   428     73        0 ( 973/1702-1009/1779) -56.3 E27-E23 N32-N38
      3  871 1693 (N28E31 L=294)   524    156        0 ( 792/1644- 928/1719)   0.0 E36-E24 N26-N29
      4 1263 1430 (N22W03 L=328)   393    139        0 (1202/1389-1315/1457)  -8.1 E01-W06 N21-N23
      5 1147 1337 (N13E01 L=323)   242    110        0 (1128/1290-1172/1392) -51.3 E03-W01 N11-N16
      6  864 1230 (N00E14 L=311)   452    156        0 ( 821/1172- 910/1295) -35.5 E17-E10 S02-N03
      7 1772 1017 (N14W48 L= 13)  1374    377        1 (1678/ 858-1828/1151) -32.3 W36-W56 N07-N20
      8 1648 1068 (N13W36 L=  1)   695    137        0 (1630/1003-1673/1133) -82.9 W34-W38 N09-N17
      9  714 1083 (S12E19 L=306)   478    222        1 ( 620/1062- 828/1112)  17.1 E25-E12 S13-S09
     10  703  913 (S21E16 L=309)   195    129        0 ( 690/ 864- 734/ 970) -81.9 E17-E14 S24-S17
     11  860  862 (S20E05 L=320)   308    112        0 ( 807/ 847- 917/ 875)  23.2 E08-E01 S21-S18
     12  625  834 (S28E19 L=305)   477    159        1 ( 548/ 817- 694/ 868)  26.6 E25-E15 S30-S24
     13 1109  747 (S20W13 L=338)   209     58        0 (1072/ 728-1139/ 763)   0.0 W10-W14 S20-S19
     14  798  641 (S35E03 L=322)  1029    240        1 ( 690/ 607- 902/ 683)  33.7 E10-W03 S38-S30
     15 1485  391 (S28W55 L= 20)  1574    173        2 (1396/ 350-1568/ 439)  48.4 W49-W59 S32-S22
     16 1252  248 (S44W52 L= 17)   652     87        1 (1222/ 226-1293/ 278)  53.1 W47-W55 S45-S41
  

    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