Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  902 1778 (N36E35 L= 76)  1547    322        6 ( 805/1700- 969/1839) -14.0 E45-E25 N31-N38
      2  567 1565 (N13E44 L= 66)   303     52        0 ( 538/1542- 602/1588)  45.0 E46-E42 N10-N14
      3 1523 1238 (N20W23 L=133)   325     93        0 (1505/1185-1536/1275) -59.0 W21-W24 N17-N22
      4 1385 1226 (N15W15 L=125)   198     58        0 (1370/1197-1405/1260) -45.0 W13-W16 N13-N16
      5  671  915 (S22E17 L= 93)   229     75        1 ( 629/ 890- 714/ 933)   9.5 E20-E14 S22-S20
      6  256  872 (S35E49 L= 61)  1583    398        2 ( 153/ 828- 365/ 911)  18.4 E60-E39 S38-S30
      7  402  545 (S53E35 L= 75)   345     60        1 ( 383/ 520- 418/ 568)  38.7 E38-E33 S55-S51
  

    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