Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1295 1784 (N41W53 L=170)   377     64        0 (1279/1754-1311/1816)  18.4 W49-W58 N39-N42
      2  241 1048 (N26E55 L= 62)  6047    819       15 ( 129/ 777- 466/1329) -85.0 E72-E36 N10-N44
      3  260  956 (N20E54 L= 63)  1040    141        3 ( 212/ 889- 285/1031)  76.0 E60-E51 N16-N24
      4 1760  928 (S21W45 L=162)   813    170        0 (1719/ 868-1802/1004)  26.6 W40-W50 S22-S17
      5 1503  897 (S14W24 L=141)   896    207        0 (1396/ 852-1588/ 946)  -8.1 W16-W30 S15-S12
      6 1732  624 (S40W40 L=157)  1532    300        1 (1620/ 543-1829/ 731)   9.5 W27-W52 S42-S35
      7  877  400 (S26E28 L= 89)  3496    915        1 ( 701/ 254-1161/ 631)  11.8 E45-E02 S35-S18
      8 1222  518 (S28E03 L=114)   680    244        0 (1149/ 447-1258/ 626)  36.0 E09-W02 S31-S22
  

    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