Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  827 1616 (N47E08 L=134)   469    100        4 ( 802/1580- 852/1650) -80.5 E10-E05 N43-N49
      2  796 1554 (N42E11 L=131)   256     49        1 ( 783/1526- 809/1583)  76.0 E12-E09 N40-N44
      3  439 1402 (N35E40 L=102)   959    275        3 ( 410/1291- 472/1511) -61.9 E45-E36 N27-N42
      4 1455 1447 (N26W36 L=178)   373     89        1 (1439/1401-1475/1489)  59.0 W33-W37 N23-N28
      5 1017  997 (N04E01 L=141)   423    139        1 ( 950/ 971-1086/1042) -29.7 E04-W03 N02-N07
      6  267  939 (N06E55 L= 87)   632    165        1 ( 227/ 918- 331/ 963)   5.2 E59-E48 N04-N07
      7  444  727 (S08E42 L=100)   374    123        0 ( 406/ 704- 479/ 758) -31.0 E44-E39 S09-S05
      8  944  496 (S27E12 L=130)   366    112        0 ( 915/ 456- 975/ 542)  45.0 E14-E09 S29-S23
  

    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