Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  959 1553 (N26E12 L=300)  1569    410        1 ( 840/1397-1040/1718) -43.6 E25-E04 N16-N36
      2 1100 1644 (N34E05 L=307)   572    128        0 (1043/1604-1154/1675) -18.4 E09-W00 N32-N36
      3 1705 1361 (N26W43 L=355)   735    122        0 (1679/1301-1727/1418) -90.0 W40-W45 N21-N29
      4 1625 1243 (N16W36 L=348)   674     94        2 (1588/1207-1666/1283) -36.9 W32-W39 N13-N18
      5  142 1029 (S14E68 L=244)   860    106        1 ( 122/ 986- 162/1086)  71.6 E71-E64 S16-S10
      6  221 1011 (S15E57 L=255)  1217    185        4 ( 151/ 957- 273/1080)  24.8 E66-E52 S18-S10
      7  488  901 (S20E33 L=279)   771    305        1 ( 411/ 837- 609/ 958)  21.8 E39-E24 S23-S14
      8  199  809 (S27E63 L=249)  1336    178        2 ( 149/ 742- 269/ 861)  18.4 E74-E53 S31-S23
      9 1178  688 (S24W15 L=327)  1764    521        0 ( 976/ 541-1338/ 841)  53.5 W04-W21 S35-S12
     10  724  806 (S23E16 L=296)   227     94        0 ( 672/ 791- 767/ 826)  26.6 E19-E13 S24-S21
     11 1505  566 (S25W41 L=353)   675    174        0 (1417/ 529-1574/ 610)  41.2 W35-W45 S29-S21
     12  577  312 (S61E37 L=275)   889     81        3 ( 542/ 284- 619/ 344)   8.7 E43-E29 S64-S58
  

    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