Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1199 1578 (N39W09 L=303)  7128   1320       11 ( 649/1397-1745/1743)  -9.1 E43-W57 N29-N48
      2  457 1368 (N19E42 L=251)   405     78        1 ( 438/1331- 483/1402) -63.4 E44-E40 N17-N21
      3  747 1195 (N11E18 L=275)  1743    441        1 ( 723/1002- 784/1389) -85.0 E22-E16 S01-N22
      4 1576 1290 (N21W37 L=330)   560    189        0 (1525/1232-1604/1373)  66.0 W32-W39 N17-N26
      5  187 1204 (N07E66 L=228)  3273    404        6 ( 150/1031- 214/1348) -81.0 E79-E60 S04-N15
      6 1499 1303 (N22W31 L=324)   326     74        0 (1480/1264-1523/1340) -59.0 W29-W32 N19-N24
      7  403 1321 (N16E46 L=248)   317     99        1 ( 354/1299- 447/1339)  14.0 E49-E41 N14-N17
      8  293 1229 (N09E54 L=240)  1113    156        3 ( 269/1162- 323/1301) -82.9 E56-E49 N05-N13
      9  345  773 (S18E47 L=247)  1952    429        1 ( 308/ 585- 376/ 966)  78.2 E52-E42 S30-S06
     10  757  647 (S23E15 L=278)   499    163        0 ( 682/ 611- 827/ 685)  26.6 E20-E10 S26-S20
     11 1526  537 (S26W40 L=334)  4701    764        9 (1165/ 394-1789/ 668)  24.1 W15-W62 S38-S16
     12  489  321 (S52E53 L=241)   680     75        1 ( 450/ 306- 522/ 343)   9.5 E63-E45 S53-S49
  

    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