Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1533 1643 (N29W51 L=269)   440     67        1 (1493/1613-1561/1664) -36.9 W48-W53 N27-N31
      2  632 1364 (N27E20 L=198)  2414    501        4 ( 447/1243- 844/1452)  13.6 E34-E04 N22-N30
      3  289 1327 (N31E50 L=168)  1758    214        3 ( 185/1294- 399/1369) -14.0 E66-E38 N27-N35
      4 1739 1310 (N06W53 L=271)   294     34        0 (1715/1296-1763/1324)   0.0 W50-W55 N05-N06
      5 1795 1042 (S11W52 L=270)   531    118        0 (1754/1009-1822/1086)  38.7 W48-W55 S12-S07
      6 1738  923 (S17W46 L=264)   276     14        0 (1720/ 904-1758/ 941)  18.4 W44-W48 S17-S15
      7 1318  831 (S15W15 L=232)   214     73        0 (1284/ 812-1352/ 848)  14.0 W12-W16 S15-S14
      8 1500  747 (S23W26 L=244)   423    104        0 (1484/ 694-1519/ 802) -82.9 W25-W27 S27-S19
      9 1602  716 (S27W34 L=252)   544    103        0 (1551/ 687-1656/ 754)  12.5 W29-W38 S28-S25
     10  442  499 (S21E51 L=166)   414    112        0 ( 392/ 476- 478/ 540)  11.3 E56-E46 S22-S19
  

    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