Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1299 1763 (N48W47 L=328)   610     94        0 (1267/1729-1345/1783)   0.0 W41-W52 N46-N49
      2 1422 1697 (N40W51 L=332)  1240    125        3 (1380/1656-1475/1746)  82.9 W45-W55 N36-N44
      3 1384 1593 (N35W40 L=320)   708    193        1 (1323/1513-1419/1636)  66.0 W36-W43 N29-N38
      4 1322 1562 (N34W33 L=314)   319    109        0 (1295/1523-1361/1610) -81.9 W32-W35 N30-N37
      5  461 1476 (N44E40 L=241)   313     50        1 ( 433/1459- 492/1490)   0.0 E43-E36 N42-N44
      6  905 1182 (N18E05 L=276)   765    318        0 ( 763/1137-1043/1227) -31.0 E13-W03 N13-N22
      7  398 1153 (N22E42 L=239)   731    212        0 ( 341/1063- 428/1218) -57.5 E47-E39 N16-N27
      8 1826  949 (S14W57 L=337)  3315    565        3 (1757/ 705-1880/1151)  65.2 W48-W66 S28-S02
      9  157 1119 (N22E70 L=210)   562     47        0 ( 138/1098- 185/1141) -12.5 E74-E65 N20-N23
     10  783 1013 (N09E15 L=266)   415    127        0 ( 728/ 977- 831/1050) -50.2 E17-E12 N06-N12
     11 1422  922 (S06W23 L=303)   660    164        0 (1385/ 837-1460/ 995)  53.1 W19-W26 S10-S02
     12  541  426 (S25E47 L=234)   335     62        0 ( 519/ 401- 559/ 457)  33.7 E49-E43 S26-S22
     13 1068  197 (S54E18 L=263)   441     64        1 (1043/ 181-1102/ 221)   0.0 E20-E13 S56-S51
  

    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