Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1759 1485 (N33W51 L= 17)   672    135        2 (1729/1454-1781/1524)  45.0 W47-W53 N30-N36
      2  295 1375 (N11E55 L=271)  2162    377        6 ( 240/1286- 366/1488) -40.9 E63-E46 N05-N18
      3 1866 1335 (N25W60 L= 26)   672     98        1 (1837/1306-1895/1369) -36.9 W56-W63 N23-N26
      4 1751 1287 (N19W47 L= 13)   483    117        1 (1708/1258-1799/1312) -14.0 W43-W51 N18-N20
      5 1450 1280 (N14W24 L=350)   726    167        3 (1379/1253-1521/1309)  12.5 W19-W28 N12-N16
      6 1592 1142 (N08W34 L=  0)   565    179        1 (1524/1113-1661/1170)   0.0 W29-W39 N06-N09
      7  443  989 (S12E37 L=289)   262     60        0 ( 411/ 974- 473/1006)   0.0 E39-E34 S12-S10
      8  332  906 (S17E46 L=280)   342     68        1 ( 298/ 893- 366/ 919)   0.0 E48-E42 S18-S16
      9 1405  865 (S11W25 L=351)   354    105        0 (1356/ 834-1464/ 895) -14.0 W21-W29 S11-S09
     10 1013  710 (S24W02 L=328)  1274    373        1 ( 854/ 615-1197/ 828)  38.3 E07-W12 S31-S15
     11  697  569 (S36E19 L=307)   467    141        2 ( 633/ 556- 774/ 585)  15.3 E23-E12 S37-S34
  

    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