Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1305 1288 (N14W23 L=279)  1450    473        3 (1172/1211-1471/1320) -14.0 W14-W34 N09-N17
      2 1837 1032 (S17W54 L=309)  1350    208        3 (1793/ 988-1897/1093)  33.7 W48-W60 S20-S13
      3  635  996 (N15E23 L=232)  1989    469        5 ( 454/ 958- 830/1041) -22.6 E36-E12 N08-N19
      4 1714  830 (S24W38 L=294)   793    171        2 (1668/ 750-1740/ 914)  37.9 W33-W42 S28-S19
      5  328  723 (N05E54 L=202)  1601    288        4 ( 290/ 619- 384/ 871)  68.2 E57-E48 S01-N14
      6 1190  808 (S09W04 L=259)   548    161        0 (1135/ 727-1250/ 861)  24.0 E01-W08 S12-S07
      7 1276  782 (S13W08 L=263)   385    145        0 (1244/ 710-1318/ 858)  40.6 W04-W12 S16-S10
      8 1685  652 (S35W34 L=289)   515    121        2 (1639/ 634-1726/ 671) -12.5 W29-W38 S37-S32
      9 1289  556 (S28W03 L=258)   337     70        0 (1255/ 526-1319/ 590)  21.8 W00-W05 S28-S26
     10  705  412 (S21E38 L=217)  2162    381        4 ( 609/ 265- 786/ 535)  79.2 E43-E35 S32-S11
     11  577  349 (S22E51 L=204)   634    213        0 ( 557/ 288- 601/ 406)  45.0 E54-E45 S27-S18
  

    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