Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  525 1705 (N24E55 L=299)  1024    200        1 ( 484/1633- 576/1799) -20.8 E74-E45 N19-N32
      2  718 1744 (N31E42 L=312)   951    153        3 ( 656/1697- 784/1795) -14.9 E49-E34 N28-N33
      3 1111 1507 (N23E06 L=347)   268     35        0 (1084/1484-1136/1531)   0.0 E07-E04 N21-N24
      4 1460 1479 (N31W17 L= 11)   519    133        0 (1404/1441-1521/1519)  -5.7 W12-W22 N29-N32
      5 1218 1447 (N22W02 L=355)   214    106        0 (1207/1395-1230/1499) -74.1 W00-W02 N18-N25
      6 1780 1184 (N21W47 L= 40)  1550    420        1 (1665/1055-1877/1356) -23.7 W34-W59 N16-N29
      7  393  963 (S23E40 L=314)   291     72        0 ( 361/ 938- 419/ 995)  45.0 E42-E38 S24-S20
      8  517  856 (S27E29 L=325)   630    191        1 ( 445/ 818- 579/ 933)  45.0 E34-E25 S30-S21
      9  647  703 (S34E17 L=336)   796    274        0 ( 543/ 651- 771/ 765)  34.5 E25-E09 S39-S27
     10  954  717 (S26W03 L=357)   192     59        0 ( 928/ 702- 979/ 734)  18.4 W01-W04 S27-S24
     11  318  563 (S51E54 L=300)   392     62        1 ( 301/ 540- 346/ 593)  33.7 E58-E49 S52-S48
     12  944  399 (S47W13 L=  7)   294     70        1 ( 906/ 387- 986/ 416)  18.4 W10-W16 S48-S44
  

    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