Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1634 1588 (N33W51 L=287)  3766    398        5 (1474/1493-1784/1662)  -9.2 W35-W67 N26-N38
      2 1597 1364 (N17W40 L=276)  1267    308        0 (1503/1250-1693/1482) -53.7 W34-W45 N10-N25
      3  407 1408 (N23E43 L=193)   306    100        0 ( 389/1361- 429/1457)  81.9 E44-E42 N19-N26
      4 1196 1392 (N19W12 L=248)   634    195        0 (1111/1340-1282/1433) -31.0 W06-W16 N15-N22
      5  325  810 (S13E49 L=187)  1136    218        1 ( 297/ 705- 354/ 918)  72.9 E52-E47 S19-S06
      6  811  784 (S17E14 L=222)   347    146        1 ( 746/ 765- 871/ 808)  -8.1 E17-E10 S17-S15
      7 1862  763 (S19W66 L=302)   868    115        1 (1830/ 735-1900/ 807) -18.4 W61-W74 S20-S16
      8 1799  614 (S28W65 L=301)   532     46        1 (1775/ 593-1823/ 634)  33.7 W61-W68 S29-S26
      9 1470  502 (S37W34 L=270)  1840    428        4 (1318/ 433-1627/ 631)  28.6 W21-W44 S41-S27
     10 1679  462 (S39W59 L=295)   705     46        1 (1658/ 430-1694/ 497) -35.5 W55-W64 S41-S36
     11  588  263 (S53E53 L=183)  1115    131        1 ( 532/ 227- 637/ 309)  11.3 E71-E41 S56-S49
  

    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