Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  967 1564 (N26E12 L=299)  1583    421        2 ( 852/1391-1057/1715) -42.1 E24-E03 N16-N36
      2 1107 1647 (N35E04 L=307)   596    190        0 (1013/1603-1171/1702) -22.6 E12-W00 N32-N37
      3 1717 1360 (N26W44 L=355)   721    155        0 (1685/1300-1739/1425) -82.9 W41-W46 N21-N29
      4 1641 1238 (N16W37 L=348)   705    104        2 (1609/1203-1682/1286) -51.3 W34-W40 N13-N19
      5  150 1031 (S14E67 L=244)   996    106        3 ( 128/ 986- 173/1088)  63.4 E70-E63 S16-S10
      6  212  987 (S16E59 L=253)   695     74        2 ( 169/ 956- 246/1022)  18.4 E64-E54 S18-S14
      7  237  842 (S25E57 L=254)  1733    358        4 ( 153/ 738- 347/ 974)  27.3 E73-E44 S31-S16
      8  498  899 (S20E32 L=279)   766    302        3 ( 424/ 835- 621/ 955)  18.4 E38-E23 S23-S15
      9 1199  692 (S23W16 L=327)  1575    525        1 ( 992/ 537-1356/ 836)  52.0 W05-W23 S36-S12
     10  741  803 (S23E15 L=296)   280    122        0 ( 672/ 787- 791/ 827)  20.6 E19-E11 S24-S20
     11 1536  572 (S24W43 L=354)   514    118        1 (1478/ 532-1582/ 605)  50.2 W40-W45 S27-S21
     12  582  310 (S61E36 L=275)   704     41        0 ( 554/ 285- 613/ 340)   0.0 E41-E31 S63-S58
  

    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