Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1732 1539 (N26W62 L=  3)  1410    176        2 (1656/1476-1800/1573) -18.4 W53-W68 N22-N28
      2 1173  829 (S15W08 L=309)   579    211        0 (1104/ 746-1235/ 916) -57.5 W04-W11 S19-S08
      3 1629  845 (S16W39 L=340)   558    163        0 (1570/ 777-1681/ 899)  41.2 W34-W43 S19-S12
      4  655  782 (S14E25 L=276)   291     88        1 ( 617/ 746- 698/ 798)  18.4 E27-E21 S15-S13
      5  566  759 (S15E31 L=270)   218     31        1 ( 544/ 738- 586/ 776) -56.3 E32-E29 S16-S13
      6  512  743 (S15E35 L=266)   385     67        2 ( 488/ 707- 542/ 769)  53.1 E37-E33 S17-S13
      7  601  624 (S23E31 L=270)   399    130        0 ( 547/ 581- 632/ 669)  32.0 E36-E28 S25-S20
      8  241  626 (S20E65 L=236)   608     53        1 ( 219/ 591- 265/ 656) -90.0 E68-E61 S22-S18
      9 1721  613 (S31W52 L=354)   444     79        0 (1681/ 577-1755/ 645)  33.7 W48-W55 S33-S29
     10  976  544 (S31E07 L=294)   288     38        0 ( 951/ 522-1001/ 562) -14.0 E08-E04 S32-S30
     11  815  288 (S49E27 L=274)   299     54        0 ( 781/ 271- 843/ 303)   8.1 E30-E23 S50-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