Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1473 1679 (N35W56 L=304)   967    167        1 (1408/1643-1550/1710) -28.6 W50-W61 N31-N37
      2  562 1651 (N59E29 L=219)   375     51        1 ( 536/1631- 585/1672)  21.8 E32-E26 N56-N60
      3 1261 1468 (N28W26 L=275)   217     56        0 (1237/1447-1292/1484)   0.0 W23-W28 N26-N28
      4  352 1343 (N39E47 L=202)   831    167        0 ( 299/1278- 389/1378) -38.7 E54-E41 N33-N41
      5  646  904 (N07E25 L=223)   233     66        0 ( 625/ 873- 668/ 933)  56.3 E27-E23 N05-N08
      6  507  879 (N08E35 L=213)   243     70        0 ( 475/ 851- 539/ 898)  18.4 E38-E32 N06-N08
      7 1744  820 (S22W45 L=294)   393     54        2 (1716/ 792-1773/ 849)   0.0 W42-W48 S23-S20
      8 1714  649 (S33W43 L=291)   365     46        2 (1690/ 628-1737/ 676)   0.0 W40-W45 S34-S30
      9  455  561 (S11E49 L=200)   594     74        1 ( 419/ 536- 497/ 588)   7.1 E52-E44 S12-S09
     10  583  442 (S21E44 L=205)   421     97        2 ( 538/ 423- 637/ 455)   0.0 E47-E38 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