Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1132 1537 (N36W02 L=292)   263    102        0 (1099/1497-1164/1581) -51.3 W00-W04 N33-N38
      2  414 1437 (N21E50 L=240)   591    133        2 ( 385/1389- 440/1493) -60.3 E54-E46 N17-N24
      3  844  860 (S10E10 L=280)   434     88        1 ( 818/ 825- 869/ 890)  63.4 E11-E07 S12-S08
      4  592  795 (S17E26 L=264)   237     51        0 ( 573/ 771- 608/ 830) -76.0 E27-E24 S18-S14
      5  881  752 (S17E06 L=284)   534    188        2 ( 798/ 730- 947/ 796)  -5.7 E12-E02 S18-S14
      6 1542  714 (S13W38 L=328)   808    190        0 (1513/ 641-1586/ 782)  90.0 W36-W41 S17-S08
      7 1308  699 (S16W22 L=312)   474    181        1 (1242/ 657-1359/ 767)  54.5 W17-W24 S19-S11
      8 1702  645 (S16W54 L=344)   447     44        1 (1680/ 622-1724/ 673)  90.0 W52-W57 S17-S14
      9  742  545 (S32E15 L=275)   238     89        1 ( 698/ 530- 772/ 569)   9.5 E18-E12 S33-S30
     10 1485  359 (S38W50 L=340)   998    110        1 (1446/ 319-1527/ 409)  90.0 W44-W53 S41-S33
  

    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