Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1399 1456 (N33W19 L= 15)  1994    438        4 (1244/1323-1592/1550) -18.4 W07-W34 N26-N37
      2  959 1482 (N27E12 L=344)   424    107        2 ( 938/1432- 984/1532) -81.9 E13-E09 N23-N30
      3  409 1480 (N19E52 L=304)   348    105        1 ( 393/1454- 433/1519) -33.7 E54-E48 N17-N21
      4  974 1278 (N15E07 L=349)   571    186        1 ( 910/1217-1054/1354) -28.6 E12-E01 N11-N18
      5  451 1274 (N07E41 L=314)   255     66        0 ( 432/1238- 472/1310)  90.0 E42-E40 N04-N09
      6 1846 1207 (N23W59 L= 55)   865    190        1 (1798/1167-1891/1239) -15.3 W53-W64 N21-N25
      7 1359 1115 (N10W19 L= 15)   648    216        6 (1298/1027-1420/1216)  65.2 W16-W22 N04-N17
  

    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