Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  732 1448 (N32E12 L= 77)   330    127        0 ( 674/1402- 781/1470)   6.3 E16-E07 N29-N33
      2  539 1404 (N33E27 L= 62)   405    109        0 ( 491/1364- 576/1440)  23.2 E30-E23 N30-N34
      3 1698 1054 (S09W44 L=132)   362     63        0 (1678/1021-1711/1092)  63.4 W41-W45 S10-S06
      4 1888 1033 (S14W64 L=152)  1194    163        4 (1854/ 986-1934/1079)  18.4 W59-W71 S16-S10
      5  220  984 (N12E57 L= 32)   368     58        0 ( 205/ 960- 239/1019) -63.4 E58-E54 N10-N14
      6 1281  895 (S11W13 L=101)   903    291        1 (1210/ 770-1344/1012)  45.0 W06-W18 S16-S04
      7  451  907 (N04E38 L= 51)   896    168        0 ( 364/ 872- 518/ 961)   9.5 E44-E31 N02-N06
      8  687  820 (S05E24 L= 65)   264     51        2 ( 656/ 802- 718/ 842) -14.0 E25-E21 S05-S03
      9  681  670 (S14E28 L= 61)   525     97        3 ( 653/ 638- 719/ 717)  31.0 E30-E24 S15-S11
  

    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