Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1170 1684 (N41W33 L=324)   417     59        0 (1137/1663-1201/1701) -63.4 W30-W36 N39-N43
      2 1029 1527 (N36W15 L=305)  1034    265        2 ( 921/1486-1161/1552) -13.2 W06-W24 N32-N38
      3 1914 1155 (S11W69 L=359)   859    113        3 (1893/1098-1936/1200)  45.0 W65-W74 S14-S08
      4 1839 1075 (S13W56 L=347)  1379    187        1 (1803/ 981-1861/1170)  59.0 W50-W61 S17-S07
      5  363 1025 (N21E43 L=247)   355     81        1 ( 332/1005- 402/1058) -38.7 E45-E40 N19-N23
      6 1345  967 (S04W17 L=307)   634    200        0 (1315/ 865-1380/1049)  52.1 W13-W20 S09-N00
      7  339  946 (N17E47 L=244)   886    250        2 ( 246/ 892- 421/ 986) -10.6 E56-E40 N12-N18
      8  323  857 (N12E50 L=241)   586    109        0 ( 298/ 812- 348/ 897) -78.7 E51-E47 N09-N14
      9 1591  617 (S32W27 L=317)  1824    311        1 (1475/ 522-1675/ 723)  22.4 W16-W34 S36-S26
     10  439  523 (S10E52 L=239)   595     81        1 ( 398/ 488- 469/ 558) -59.0 E53-E49 S12-S07
     11  587  298 (S28E55 L=236)  1681    238        4 ( 487/ 257- 674/ 341) -23.2 E60-E44 S31-S22
     12  904  211 (S42E37 L=253)  2349    305        3 ( 798/ 126-1029/ 297)  11.6 E61-E22 S48-S37
  

    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