Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1604 1650 (N22W64 L=209)  1176    116        1 (1560/1606-1655/1692)   0.0 W56-W76 N19-N22
      2  366 1537 (N56E45 L=101)   477     66        1 ( 340/1510- 403/1554)  -5.2 E50-E38 N54-N58
      3  800 1460 (N37E01 L=144)   326    101        2 ( 743/1438- 848/1474) -26.6 E04-W02 N34-N38
      4 1832 1401 (N01W70 L=215)  1276    136        2 (1806/1340-1856/1467)  54.5 W64-W75 S02-N05
      5  162  870 (N17E68 L= 77)  1173    116        1 ( 134/ 819- 189/ 938) -90.0 E72-E64 N13-N20
      6 1362  780 (S17W13 L=158)   849    277        0 (1259/ 717-1469/ 871)  10.6 W05-W21 S18-S14
      7 1010  750 (S09E08 L=137)   743     97        2 ( 938/ 636-1070/ 850) -86.2 E09-E06 S16-S01
      8 1442  706 (S24W16 L=161)   312    205        0 (1399/ 664-1476/ 747)  31.0 W12-W18 S26-S21
      9 1475  632 (S30W17 L=162)   420    128        1 (1422/ 604-1528/ 670)  -6.3 W12-W21 S31-S27
     10 1117  418 (S32E13 L=132)   974    197        2 (1068/ 332-1153/ 489)  71.6 E17-E10 S38-S26
  

    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