Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  730 1737 (N62W04 L=244)   533    107        1 ( 697/1703- 778/1774)  14.0 E01-W11 N58-N64
      2 1383 1696 (N34W51 L=291)   764     97        2 (1339/1661-1423/1722) -59.0 W46-W54 N31-N36
      3 1094 1633 (N40W26 L=266)   336     73        0 (1073/1604-1129/1667)  78.7 W24-W27 N37-N42
      4  509 1638 (N61E29 L=211)   359     58        0 ( 488/1612- 532/1660)  36.9 E31-E24 N59-N63
      5 1136 1472 (N29W20 L=260)   217     68        1 (1104/1456-1177/1484) -11.3 W18-W23 N28-N30
      6  440 1194 (N31E35 L=205)   373    100        0 ( 392/1175- 487/1212)  -6.3 E39-E30 N29-N31
      7 1832  915 (S23W53 L=293)   783     93        1 (1810/ 869-1853/ 956)  33.7 W50-W56 S25-S21
      8 1799  789 (S30W48 L=288)   660    126        1 (1759/ 764-1854/ 816)  -4.4 W43-W56 S32-S27
      9 1105  651 (S16E05 L=235)   320     86        2 (1064/ 632-1155/ 671) -21.8 E07-E02 S17-S14
     10 1057  498 (S25E13 L=227)   633    269        0 ( 992/ 451-1112/ 561) -29.7 E16-E09 S27-S20
     11  732  441 (S20E35 L=205)   425    121        2 ( 684/ 422- 797/ 463) -39.8 E38-E31 S22-S17
     12  630  378 (S22E46 L=194)  1651    212        3 ( 565/ 329- 721/ 432) -31.0 E52-E38 S25-S17
  

    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