Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1662 1454 (N24W47 L=243)   306     37        0 (1638/1434-1685/1473) -18.4 W44-W48 N22-N25
      2  285 1290 (N14E53 L=143)  1632    338        2 ( 256/1160- 315/1436) -74.5 E57-E49 N06-N24
      3  664 1313 (N15E23 L=173)   659    184        3 ( 593/1289- 759/1336)  -5.2 E27-E16 N13-N16
      4  217 1221 (N10E59 L=137)   517     81        0 ( 195/1182- 240/1260)  53.1 E60-E56 N08-N12
      5  328 1168 (N07E47 L=149)   852    146        2 ( 289/1104- 371/1226)  45.0 E50-E43 N02-N10
      6  884  837 (S14E09 L=187)   937    319        0 ( 820/ 697- 919/ 986)  70.6 E13-E06 S22-S05
      7 1358  794 (S17W21 L=217)   387    208        1 (1314/ 708-1392/ 851)  60.9 W18-W23 S22-S13
      8 1438  793 (S17W27 L=223)   274    108        1 (1389/ 774-1483/ 823) -18.4 W23-W29 S18-S15
      9  892  473 (S38E10 L=185)   737    212        0 ( 780/ 453- 982/ 506)   7.1 E19-E03 S39-S35
     10 1685  443 (S39W62 L=258)   571     72        2 (1648/ 424-1723/ 461)   0.0 W56-W70 S40-S38
     11  898  312 (S51E13 L=183)   296     57        1 ( 869/ 294- 921/ 336)  11.3 E15-E10 S52-S48
  

    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