Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1188 1735 (N53W23 L= 51)   350     71        0 (1150/1721-1227/1748)   6.3 W18-W27 N51-N53
      2 1575 1565 (N35W51 L= 79)  1552    309        2 (1428/1527-1697/1632) -23.5 W38-W64 N30-N41
      3  556 1426 (N32E33 L=355)   927    219        2 ( 498/1366- 621/1519) -33.7 E41-E26 N28-N38
      4 1371 1302 (N20W25 L= 53)   223    100        0 (1351/1270-1395/1350)  78.7 W23-W26 N17-N22
      5  558 1204 (N18E30 L=358)   382    155        0 ( 522/1140- 612/1249)  54.5 E32-E26 N13-N20
      6  705 1150 (N14E20 L=  8)   530    249        1 ( 615/1067- 770/1225) -42.3 E26-E15 N08-N18
      7 1066 1105 (N09W03 L= 31)   450     90        1 (1029/1066-1094/1145) -59.0 W00-W04 N06-N12
      8 1410  786 (S13W24 L= 52)   487    118        1 (1363/ 753-1467/ 821)  15.9 W20-W27 S14-S10
      9  430  750 (S10E42 L=346)   869    152        0 ( 379/ 705- 500/ 780)  20.0 E47-E36 S13-S08
  

    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