Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  975 1358 (N17E10 L=283)  3311    788        8 ( 749/1068-1184/1645)  15.7 E32-W08 N03-N32
      2 1630 1586 (N49W39 L=331)   937    168        1 (1579/1549-1676/1609)  10.3 W33-W45 N46-N51
      3 1812 1072 (N17W54 L=346)  2037    366        3 (1721/ 927-1907/1193) -25.6 W44-W69 N10-N23
      4  261 1065 (S13E53 L=239)   328    102        2 ( 243/1021- 280/1107) -68.2 E55-E51 S15-S10
      5 1203 1001 (N01W11 L=303)   249     63        1 (1178/ 977-1231/1039) -33.7 W09-W12 N00-N02
      6  806  842 (S16E10 L=282)  1561    443        2 ( 607/ 791-1016/ 904)  22.6 E22-W02 S21-S08
      7  350  861 (S24E43 L=249)   958    365        2 ( 239/ 832- 488/ 897)  10.3 E53-E31 S25-S21
      8 1231  405 (S35W30 L=322)  4514    578       10 ( 868/ 365-1505/ 479)  24.9 W03-W46 S45-S24
  

    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