Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1203 1780 (N46E16 L=330)  1479    173        5 (1158/1693-1247/1841) -32.7 E22-E07 N40-N51
      2  703 1645 (N21E39 L=306)   400     56        0 ( 673/1624- 734/1668)   0.0 E42-E36 N19-N21
      3 1064 1577 (N26E14 L=331)   299     76        1 (1024/1553-1103/1599)  51.3 E16-E12 N23-N28
      4 1668 1509 (N44W32 L= 17)  1471    234        3 (1590/1426-1752/1594) -18.4 W23-W41 N39-N48
      5 1111 1470 (N21E08 L=338)   357    134        1 (1096/1401-1123/1526) -63.4 E10-E05 N16-N24
      6 1513 1403 (N29W19 L=  5)   870    264        0 (1466/1268-1535/1499) -55.3 W13-W22 N21-N34
      7 1279 1367 (N19W05 L=350)   327    131        2 (1227/1324-1302/1401) -18.4 W01-W07 N17-N21
      8  348 1030 (S23E44 L=301)   405    130        1 ( 331/ 988- 381/1080) -59.0 E46-E41 S25-S20
      9  737 1016 (S15E16 L=329)   228     75        1 ( 715/ 986- 775/1039)  63.4 E17-E14 S16-S12
     10  583  885 (S26E23 L=323)   392    160        1 ( 510/ 854- 653/ 918)  37.9 E27-E18 S29-S22
     11 1096  651 (S25W16 L=  1)   776    285        0 ( 992/ 586-1231/ 692)  47.5 W11-W22 S31-S19
  

    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