Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  757 1710 (N58E08 L=  1)  1841    266        4 ( 674/1640- 884/1756)  -4.6 E18-W07 N53-N61
      2 1204 1701 (N46W33 L= 42)  1109    181        0 (1121/1672-1292/1724)  -6.7 W25-W42 N43-N48
      3 1439 1675 (N38W51 L= 60)   825    120        1 (1390/1650-1501/1690) -15.3 W46-W58 N35-N39
      4 1598 1441 (N20W49 L= 58)   579    105        3 (1552/1394-1643/1468) -33.7 W45-W52 N16-N21
      5  597 1338 (N33E24 L=345)   300     74        0 ( 568/1310- 636/1354)   0.0 E26-E21 N31-N33
      6 1320 1262 (N15W23 L= 32)  1164    262        3 (1248/1149-1363/1354) -90.0 W20-W24 N08-N21
      7  751 1218 (N23E14 L=355)   392    134        1 ( 734/1151- 764/1284)  83.7 E15-E12 N18-N27
      8  584 1080 (N17E27 L=342)   497    156        0 ( 513/1052- 648/1111)   5.7 E32-E22 N15-N18
      9 1877 1084 (S08W65 L= 74)   526     53        1 (1858/1055-1898/1109)  56.3 W62-W68 S09-S06
     10 1741  711 (S28W47 L= 56)  2456    233        6 (1698/ 605-1774/ 814) -69.4 W42-W52 S35-S19
     11  928  220 (S47E30 L=339)   511     27        1 ( 903/ 203- 956/ 242) -26.6 E32-E27 S49-S44
  

    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