Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  199 1274 (N20E68 L=200)  1105    103        6 ( 168/1235- 237/1319) -45.0 E74-E62 N17-N23
      2 1647 1232 (N14W44 L=312)   901    225        1 (1560/1163-1718/1305) -39.3 W38-W49 N09-N18
      3  420  837 (S06E41 L=226)   300     84        1 ( 377/ 824- 453/ 856)   8.1 E45-E38 S06-S04
      4  542  703 (S15E34 L=234)   238     47        0 ( 523/ 676- 565/ 730)  45.0 E35-E31 S16-S12
      5  418  592 (S22E47 L=221)   671    204        0 ( 378/ 507- 456/ 663)  74.7 E52-E44 S28-S17
      6 1240  625 (S22W13 L=281)   305     64        0 (1220/ 596-1269/ 648)  45.0 W11-W14 S23-S20
      7 1142  529 (S28W06 L=274)   253     44        0 (1123/ 506-1163/ 559)  90.0 W04-W07 S29-S25
      8 1440  390 (S41W32 L=300)   436     45        0 (1412/ 370-1465/ 412)  26.6 W30-W34 S42-S39
  

    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