Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1444 1646 (N47W10 L=356)   616    127        2 (1396/1621-1486/1678)  40.6 W06-W13 N44-N51
      2  731 1579 (N20E34 L=312)   569    184        1 ( 661/1548- 791/1626)  18.4 E39-E30 N17-N23
      3 1148 1301 (N14E00 L=345)  1213    355        0 (1069/1128-1255/1457) -37.9 E09-W09 N07-N21
      4  799 1342 (N07E21 L=325)   296    122        0 ( 759/1311- 833/1373)  45.0 E23-E19 N04-N09
      5 1702 1280 (N30W37 L= 23)   500    116        1 (1665/1215-1725/1319) -45.0 W33-W39 N26-N32
      6  205 1207 (S14E62 L=284)  1123    212        1 ( 152/1167- 267/1244)  20.0 E68-E55 S16-S11
      7 1810 1045 (N19W51 L= 37)   393     74        1 (1792/1008-1831/1078) -36.9 W48-W53 N17-N21
      8 1904  986 (N20W63 L= 49)   580    123        2 (1880/ 952-1925/1023)  63.4 W60-W65 N16-N22
      9 1728  921 (N10W47 L= 33)   342     42        1 (1710/ 895-1744/ 950)  90.0 W45-W48 N07-N11
     10 1253  756 (S13W20 L=  6)   335    153        0 (1170/ 724-1324/ 801)  52.1 W16-W23 S17-S08
     11 1844  741 (N04W67 L= 53)   742    142        4 (1814/ 701-1868/ 776) -26.6 W62-W71 N02-N05
     12  456  691 (S40E30 L=316)   521    149        0 ( 393/ 657- 517/ 733)  -5.2 E35-E24 S40-S38
     13 1511  613 (S13W42 L= 28)   260     54        0 (1480/ 596-1532/ 637)  56.3 W40-W43 S14-S11
     14 1710  584 (S09W60 L= 46)   920    172        0 (1636/ 547-1755/ 636)  48.8 W54-W63 S12-S04
     15 1454  432 (S25W48 L= 33)  1262    228        0 (1347/ 405-1543/ 463)  23.2 W39-W53 S28-S20
     16 1334  277 (S38W50 L= 36)   425     94        2 (1308/ 259-1359/ 309) -18.4 W47-W53 S39-S36
     17 1214  273 (S43W43 L= 29)   466     71        1 (1182/ 258-1245/ 293)   0.0 W39-W45 S44-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